
Calhoun 

iniutuiiaiul AKliiv« ou tfit Nilvdl Poi($ra{jua(« School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 



Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


1998-09 

Development and verification of an 
aerodynamic model for the NPS FROG UAV 
using the CMARC panel code software suite 

Pollard, Stephen J. 

Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/8264 


Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


http://www.nps.edu/ljbrary 


CsMwun is the Neval Postgraduate School's public access distal repository for 
research oiateriels and tnstitutjiooal pubftcatiions created by the NPS community. 
Cathouni is named for Professor of Mathematcs Guy K. CatHiuo, NPS's first 
appointed — and publi^d — scholar^ author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
MontereVr California USA 93943 



NPS ARCHIVE 
1998.09 
POLLARD, S. 


NAVAL POSTGRADUATE SCHOOL 
Monterey, California 



THESIS 


Thesis 

P 6875 


DEVELOPMENT AND VERIFICATION OF AN 
AERODYNAMIC MODEL FOR THE NPS FROG UAV 
USING THE CMARC PANEL CODE SOFTWARE SUITE 

by 

Stephen J. Pollard 
September 1998 

Thesis Advisor; Max F. Platzer 

Co-Advisor/Second Reader: Kevin D. Jones 


Approved for public release; distribution is unlimited. 




DUDLEY KNOX LlBRA.Rf 
mavAL ^C$TG1VsD»JA1 F SCHCOi 
'■.-J'NT rjRE^ CA A''\Cl 


REPORT DOCUMENTATION PAGE 

Form Approved 

0MB No. 0704-01c.> 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instruction 

searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Sem 
comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, 
Washington headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, v/ 
22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188) Washington DC 20503. 

1. AGENCY USE ONLY fLeave bfanfr; 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

September 1998 Aeronautical Engineer Thesis 

4. TITLE AND SUBTITLE 

DEVELOPMENT AND VERIFICATION OF AN AERODYNAMIC MODEL FOR THE NPS 
FROG UAV USING THE CMARC PANEL CODE SOFTWARE SUITE 

5. FUNDING NUMBERS 

6. AUTHOR(S) 

Pollard, Stephen J. 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Naval Postgraduate School 

Monterey, CA 93943-5000 

8. PERFORMING 
ORGANIZATION REPORT 
NUMBER 

9. SPONSORING / MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

10. SPONSORING/ 
MONITORING 

AGENCY REPORT NUMBt R 


11. SUPPLEMENTARY NOTES 


The views expressed in this thesis are those of the author and do not reflect the official policy or position of the Department of 
Defense or the U.S. Government. 


12a. DISTRIBUTION / AVAILABILITY STATEMENT 


12b. DISTRIBUTION CODE 


!! 

!l 


Approved for public release; distribution is unlimited. I 

13. ABSTRACT (maximum 200 words) j 

The CMARC panel-code is evaluated for the development of an aerodynamic model of the Naval Postgraduat e 
School FROG Unmaimed Air Vehicle (UAV). CMARC is a personal computer hosted panel-code software suite for 
solving inviscid, incompressible flow over complex three-dimensional bodies. A panel model of the NPS FROG UAV I 
is developed to obtain stability derivative data at the cruise flight condition. Emphasis is placed on comparing the [ 
CMARC data to aerodynamic models obtained from classical design techniques and parameter estimation. Linearized ' 
longitudinal and lateral-directional state-equation models are used to compare the dynamic response of each data set. 
In addition, CMARC is used to generate static-source and angle-of-attack sensor position corrections. Position 
corrections are provided in look-up table and curve-fit formats. The aerodynamic model obtained with CMARC I 
demonstrated higher fidelity dynamic longitudinal response than the classical design model. Dynamic lateral- \ 
directional response is similar to that obtained from classical design techniques. Adjustment through comparison with I 
flight-test data is still required to optimize the CMARC model. Future studies should concentrate on improving i 
CMARC modeling of fuselage side force through the addition of wake separation lines. Additionally, the propellci 
disk should be modeled in an attempt to capture the effects of increased dynamic pressure over the horizontal ard 


vertical tail surfaces. 


14. SUBJECT TERMS 

Unmanned Aerial Vehicles, UAV, CMARC, Panel Method Ames Research Center, PMARC, Panel Code, 
Stability Derivatives, Boundary Layer Code, Aircraft D 3 mamic Response 

15. NUMBER 0? 
PAGES 

195 

16. PRICE CODE 

17. SECURITY CLASSIFICATION 

OF REPORT 

Unclassified 

18. SECURITY CLASSIFICATION OF 
THIS PAGE 

Unclassified 

19. SECURITY CLASSIFICATION 

OF ABSTRACT 

Unclassified 

20. LIMITATION 
OF ABSTRACT 

UL 


NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 


Prescribed by ANSI ft i 


1 


































' ft - ' , ^ f 

V ' .. ' 

■ , - '^SSOF, 

% " . ^ - ! *■*■ ¥1 . “Hll . 

. iriJ 

JIL ( 11 

■ .Jli- 

■■ j 

<S.,- ^ 

'K 

i 


rt’. 



11 






Approved for public release; distribution is unlimited 


DEVELOPMENT AND VERIFICATION OF AN AERODYNAMIC MODEL FOR 
THE NPS FROG UAV USING THE CMARC PANEL CODE SOFTWARE SUITE 

Stephen J. Pollard 
Commander, United States Navy 
M.S., Naval Postgraduate School, 1997 

Submitted in partial fulfillment of the 
requirements for the degree of 


AERONAUTICAL AND ASTRONAUTICAL ENGINEER 


from the 


NAVAL POSTGRADUATE SCHOOL 
September 1998 



OUDLF.Y KNOX LiBRABY 
^!•AVAL POSTGRADiJATE SChC' 

JviOvTFi^HY CA 9:/>3-A101 

ABSTRACT 


The CMARC panel-code is evaluated for the development of an aerodynamic 
model of the Naval Postgraduate School FROG Unmanned Air Vehicle (UAV). 
CMARC is a personal computer hosted panel-code software suite for solving inviscid, 
incompressible flow over complex three-dimensional bodies. A panel model of the NPS 
FROG UAV is developed to obtain stability derivative data at the cruise flight condition. 
Emphasis is placed on comparing the CMARC data to aerodynamic models obtained 
from classical design techniques and parameter estimation. Linearized longitudinal and 
lateral-directional state-equation models are used to compare the dynamic response of 
each data set. In addition, CMARC is used to generate static-source and angle-of-attack 
sensor position corrections. Position corrections are provided in look-up table and curve- 
fit formats. The aerodynamic model obtained with CMARC demonstrated higher fidelity 
dynamic longitudinal response than the classical design model. Dynamic lateral- 
directional response is similar to that obtained from classical design techniques. 
Adjustment through comparison with flight-test data is still required to optimize the 
CMARC model. Future studies should concentrate on improving CMARC modeling of 
fuselage side force through the addition of wake separation lines. Additionally, the 
propeller disk should be modeled in an attempt to capture the effects of increased 
dynamic pressure over the horizontal and vertical tail surfaces. 
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I. INTRODUCTION 


A. BACKGROUND 

Computational fluid dynamics (CFD) is increasingly used as a design and analysis 
tool. As the price of computer hardware drops and computational power increases, CFD 
becomes more attractive to a larger audience. CFD tools range from the high-end three- 
dimensional (3D) Navier-Stokes solvers for compressible, viscous fluids to potential flow 
solvers for incompressible, inviscid flows. This paper discusses the development of an 
aerodynamic model of the Naval Postgraduate School (NPS) Fiber Optic Guided (FROG) 
Unmanned Air Vehicle (UAV) using a personal computer hosted panel code. 

Validation of the Personal Simulation Works panel code software and initial 
FROG UAV modeling was completed in June 1997. Results are reported in Ref [1], 
This study expands on the modeling effort to include the development of a complete 
aerodynamic model at the trim flight condition. Stability derivative data are obtained 
from the panel code and then compared to data from classical design calculations and 
parameter estimation. A linearized state space simulation is used to extract the air vehicle 
dynamic modes and to model control response. Results are compared for all three data 
sets. 

The Personal Simulation Works software suite, consisting of LOFTSMAN, 
CMARC and POSTMARC, is used for all aspects of the study. The software provides 
for panel model development, input file processing and the visualization of results. 
Emphasis is placed on verifying both the accuracy and suitability of the CFD programs 
for aerodynamic modeling. 

Until recently, personal computers (PC) did not have the computational power or 
memory to be practical for panel code CFD programs. Things have changed with the 
introduction of the Pentium class PC and low cost RAM. AeroLogic capitalized on the 
power of the Pentium class PC and developed Personal Simulation Works (PSW). PSW 
is based on the 3D, low order, inviscid potential flow solver named CMARC. CMARC is 
a re-hosted version of NASA’s Panel Method Ames Research Center (PMARC). 
PMARC was re-written in the C language and compiled for IBM compatible PCs. 
CMARC runs under the DOS operating system. It will also run in a DOS window under 
the WINDOWS 3.x, 95 or NT operating systems. Enhanced capabilities include; 
improved memory management, an expanded set of command line switches and 
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provisions for expanded boundary layer post-processing capabilities. However, the core 
processing algorithms remain the same as implemented in PMARC. 

LOFTSMAN, the PSW pre-processing program, is used to mesh complex 3D 
bodies and create input file patches. The program runs under the Windows operating 
system and allows the user to loft conics based 3D surfaces. The program automatically 
creates CM ARC, PMARC or VS AERO input patches based on desired panel densities 
and distribution. 

POSTMARC is used for flow visualization and integration of resultant forces. It 
runs under the Windows 3.x, 95 and NT operating systems. POSTMARC reads CMARC 
or PMARC output files and provides for the visualization of model geometry, wake 
stepping, on and off-body streamlines and surface phenomena. 

B. REQUIREMENTS 

The Naval Postgraduate School Aeronautics Department is integrating UAV 
hardware and software to demonstrate autonomous flight, trajectory tracking and 
automatic landing. Closed-loop flight control development requires a valid aerodynamic 
truth model for the UAV airframe. The introduction of each new airfi-ame requires the 
development of a new aerodynamic truth model. Most recently, Papageorgiou [Ref 2] 
developed and tested an aerodynamic model for the NPS FROG UAV based on classical 
methods. His model produced a close match to flight test results in the longitudinal axis. 
However, the lateral-directional axis required modifications based on measured flight test 
data to produce acceptable results. With the availability of low cost panel code CFD 
capabilities, it is suggested that a panel code model of the FROG UAV will provide more 
refined stability derivative data for an initial aerodynamic model. 

In addition to a valid truth model, accurate pitot-static and angle-of-attack sensors 
are required for highly augmented flight control systems. CMARC is well suited for 
solving on-body static pressure distributions and off-body flow velocities over the 
predominately attached flow fields of a fuselage forebody. This proves particularly 
useful for generating pitot-static and angle of attack correction curves and look-up tables. 

C. STATEMENT OF OBJECTIVES 

The Naval Postgraduate School Aeronautics Department has both active CFD 
research and avionics development programs. The primary purpose of this investigation 
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is to verify the accuracy and suitability of the PSW software suite for developing an 
aerodynamic model of the NPS FROG UAV. Specific objectives are as follows: 

• Demonstrate panel code modeling, processing and visualization on a Pentium 
PC using the PSW software package. 

• Develop and analyze a panel code model for the NPS FROG UAV using PSW 
to obtain a complete set of stability derivatives at the cruise trim condition. 

• Verify the aerodynamic model through comparison to data obtained from 
classical design calculations and parameter estimation. 

• Demonstrate techniques for producing angle-of-attack vane and pitot-static 
correction curves. 

• Develop a user guide that outlines PSW panel code modeling for the 
extraction of stability derivatives. 
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11. OVERVIEW OF PERSONAL SIMULATION WORKS 


A. GENERAL 

Personal Simulation Works is a PC based software suite that provides for the three 
primary CFD requirements; 3D modeling of an aircraft (LOFTSMAN), panel code flow 
solver (CMARC), and post-processing of the computed flow field (POSTMARC). The 
software package contains three applications hosted on the IBM compatible personal 
computer. Each software program is discussed separately. 

B. LOFTSMAN 

LOFTSMAN is a Windows hosted 3D modeling tool that generates surface panel 
distributions for CMARC or PMARC input files. The program is based on conics, which 
allows rapid lofting of streamlined bodies such as aircraft fuselages and engine nacelles. 
In addition, wing and control surfaces can be designed with the extensive library of airfoil 
templates or with user specified coordinates. The software is well documented, including 
a tutorial, in the Personal Simulation Works User Guide [Ref. 3]. LOFTSMAN is 
primarily designed for creating new objects, but an existing airframe can be matched 
quite closely with just a detailed three-view drawing that includes frame cross sections. 
Appendix A outlines the development of the FROG UAV model. 

1. Streamlined Bodies 

LOFTSMAN functionality is divided into Body Objects and Wing Objects. In 
general, they remain separate unless the intersection between a wing and body is required. 

Body Objects are created using a family of curves called second-degree conics. 
Circles, ellipses, parabolas and hyperbolas are among this group. An entire fuselage is 
described by specifying just four lines. These are the top waterline (TW), bottom 
waterline (BW), the maximum breadth line (MB) and the waterline of the maximum 
breadth line (WW). For each line, the beginning, ending and a few points along the line 
are specified. Control points are also specified with a curvature factor that allows 
LOFTSMAN to generate a smooth conic between the points. The power of conic lofting 
will become evident when discussing the modeling of the complex FROG UAV fuselage 
in Chapter V. 


5 


•> 


Wings and Control Surfaces 


Wings and control surfaces are easily specified in LOFTSMAN using a short 
input file created with any text editor. The file specifies root, intermediate and tip rib 
section, location, axis, chord and incidence. LOFTSMAN then fairs a smooth surface 
through the rib sections. Washout is specified by varying the incidence of the root and tip 
ribs. Sweep-back is controlled by staggering the tip rib location with respect to the root 
rib. Once the general wing surface is specified, control surfaces such as ailerons, flaps 
and elevators can be deflected and meshed. 

3. Patches 

LOFTSMAN automatically meshes 3D surfaces and creates patches for 
CM ARC/PM ARC input files. The distinction between a mesh and a patch is important. 
A mesh is a set of quadrilateral and triangular panels that represent the surface of a wing 
or body. When the set of panels is organized and formatted to create a sub-component 
portion of a CMARC or PMARC input file, it is called a patch. 

A body or wing surface is first meshed at a density specified by the user. Panel 
compression options include cosine and half-cosine spacing. After meshing the object, 
one saves it to a text file as a formatted patch. One then opens the patch file with any text 
editor and copies/pastes the patch text into the appropriate location in the CMARC input 
file. 

Each control surface deflection requires a separate mesh and formatted patch. For 
instance, to evaluate roll performance one needs to separately mesh an upward aileron 
deflection on the right wing and a downward deflection on the left wing. If multiple 
deflections of a single control surface are required, each deflection must be meshed 
separately. 

C. CMARC 

CMARC is the C version of PMARC low-order, 3D panel code. Inviscid, 
irrotational, incompressible, potential flow is assumed. Low-order means that source and 
doublet strength distribution is constant across a panel. There is no attempt to match the 
source or doublet strength of an adjacent panel at a common edge. Advanced features 
include internal flow modeling and time stepping wake models. 
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PM ARC version 12.19 was released as FORTRAN 77 source code in 1992. 
CMARC was rewritten in the C language and compiled for hosting on IBM compatible 
personal computers by AeroLogic, Inc. The program runs under the DOS operating 
system. It will also run in a DOS window from Windows 3.1, 95 or NT. Enhanced 
features include command-line options and flexible memory management. Command 
line options simplify batch processing by adding an extensive set of switches that can be 
set external to the CMARC input file. Flexible memory management provides for the 
automatic sizing of arrays without having to recompile the source code. 

D. POSTMARC 

POSTMARC is a Windows post-processing program for the visualization of 
CMARC and PMARC output files. Capabilities include body geometry, wake stepping, 
surface pressure and streamline visualization. POSTMARC also provides the capability 
to integrate pressure and skin friction forces over the model geometry. This proves 
particularly useful when one desires to recalculate loads around a different center of 
gravity. 

An interesting feature for design work is the integration of panel surface area to 
obtain total wetted area. After lofting a new geometry in LOFTSMAN, a quick check of 
geometry is made by running CMARC with the -g command line toggle. The total 
wetted area is then checked in POSTMARC. This function is particularly useful when 
working to reduce skin friction drag. 

Versions 1.17.3 and later of POSTMARC include the capability to integrate skin 
friction drag coefficient over the model geometry. It is important to note that a key piece 
of the drag equation is missing from a POSTMARC solution. CMARC provides induced 
drag from the surface pressure distribution and skin friction drag from the 2D boundary 
layer code. Skin friction is only calculated up to the point of boundary layer separation. 
Pressure drag due to separation, a major portion of the drag equation, is missing from a 
CMARC/POSTMARC solution. 

In fact, if one isn’t careful, POSTMARC drag calculations can be misleading. 
Take for instance two similar model configurations with only minor geometry differences 
that do not affect wetted area. It is possible for the model with more flow separation to 
have less skin friction drag because there is no CMARC output for skin friction 
coefficient after the boundary layer code predicts separation. During iterative design 
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work, this could lead to the incorrect conclusion that the design team is reducing overall 
drag. Perhaps a better function for LOFTSMAN than integrated skin friction drag would 
be a function that predicts the percentage of attached flow and laminar flow. Iterative 
design changes could be made that maximize laminar flow and minimize separated flow. 
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III. CMARC PANEL CODE THEORY 


A. POTENTIAL FLOW PANEL CODE THEORY (CMARC/PMARC) 

Potential flow theory involves the superposition of sources and doublets to 
generate the desired flow field around a 3D body. It assumes inviscid, irrotational and 
incompressible flow. As such, valid solutions are only obtained at low Mach numbers 
and for flow fields without large areas of separation. 

The basic concept of panel methods, as outlined by Bertin and Smith [Ref. 4], 
requires the modeling of the desired 3D configuration with a large number of 
quadrilateral and triangular panels representing the surface of the aircraft. A series of 
sources, doublets and vortices is then distributed on each panel. Superposition allows the 
simultaneous computation of the singularity strengths required to satisfy flow tangency 
on the surface. The inviscid, irrotational and incompressible flow field represented by the 
superposition of sources and doublets satisfies the Laplace equation: 


V^O = 0 


3.1 


Using Green’s Theorem, the potential at any point P in the flow is represented by: 




3.2 


Where (O-O,.) represents the potential fi'om the doublet distribution and 
n . (VO - VO.) represents the potential fi'om the source distributions. 

CMARC is a low order panel code that assumes constant source and doublet 
strength distributions across each panel. Figure 3.1 shows a panel layout for a generic 3D 
wing fuselage configuration. It is important to note that for a 3D solution, there is an 
equivalence to surface doublet and surface vortex distributions. CMARC implements 
source and doublet distributions. 
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Figure 3.1 Typical Wing/Body Panel Code Configuration, from Ref. [5]. 

As mentioned previously, the general boundary condition imposed is tangential 
flow at the surface. CMARC, as outlined in Ref [3], allows the modification of this 
boundary condition on individual panels or groups of panels. A normal surface velocity 
distribution may be specified to simulate flow into or out of ducts. 

In order to produce lift, a potential flow panel code requires a method to 
implement the Kutta condition. As noted in Anderson [Ref 6], the Kutta condition at the 
trailing edge implies that the circulation, F, around an airfoil is such that the flow exits 
the trailing edge smoothly. In addition, the velocities leaving the top and bottom surfaces 
are finite and equal in magnitude and direction. 

Panel codes impose the Kutta condition by the shedding of wake panels along the 
trailing edges or separation lines. Wake panels are similar to a surface panel with only a 
doublet distribution. The doublet strength of the attached wake panel equals the 
difference in doublet strengths of the two adjacent surface panels. 
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The CMARC core panel code processing engine is functionally equivalent to the 
PMARC panel code module. The implemented equations are well documented by Ashby 
et al. [Ref. 5]. The PMARC documentation includes a wing-body combination, shown in 
Figme 3.1, evaluated by PMARC with good correlation to experimental data. The 
results are shown in Figures 3.2 and 3.3. In addition, Lambert [Ref 7] compared 
PMARC panel code results to several theoretical and experimental test cases with good 
correlation at low angle-of-attack. Sensitivity to wake placement is highlighted by his 
studies. 

Wake positioning can have a large influence on potential flow solutions. A wake 
is obviously attached to the trailing edge of wings and control surfaces with sharp, thin 
trailing edges to produce the Kutta condition. However, wake positioning on streamlined 
fuselages, missile airframes and nacelles is more of an art than science. Recently, Tuncer 
and Platzer [Ref 8] investigated generalized wake placement techniques for cylindrical 
bodies of revolution with good correlation to experimental data at up to 20 degrees angle- 
of-attack. The techniques were used with success by this author in Ref [1] for the 
verification of CMARC calculations for flow over an inclined 6:1 prolate spheroid. 



Figure 3.2 Comparison of Experimental Data and PMARC Results for Two 
Spanwise Stations of the Wing/Body (a = 4°), from Ref. [5]. 


11 






















□ 

o exp. (ref. 13) 



upper surface 



• exp. (ref. 13) 

Cp 0 - 

X-- 

lower surface 



PMARC. upper surface 


/ • 

*•* PMARC, lower surface 

.2 ^ 

1 ■ I- ^ t ' 1 -—1-- - 1 --- \ --- \ ■ ■ 


0 

10 20 30 40 50 60 70 80 90 100 


Percent fuselage length 



Figure 3.3 Comparison of Experimental Data and PMARC Along the Fuselage 
Centerline of the Wing/Body Configuration (a = 4°), from Ref. [5]. 


B. CMARC BOUNDARY LAYER ANALYSIS THEORY 

CMARC and PMARC use the same two-dimensional integral method to calculate 
boundary layer characteristics along a surface streamline. A transition model 
automatically switches from laminar to turbulent calculations. The developers of the 
PMARC code chose a 2D integral routine over a 3D finite difference grid method due to 
speed and robustness of the calculations [Ref 5]. Building a finite difference grid is a 
difficult and time consuming process requiring the user to develop grids over complex 3D 
surfaces. In addition, boundary layer calculation times can easily exceed that required for 
the basic potential flow solution. Reference [9] gives a good outline of three-dimensional 
finite difference methods. 
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The CMARC and PMARC User’s Guides [Refs. 3 and 5] contain detailed 
discussions on the development of the CMARC/PMARC boundary layer code starting 
from the two-dimensional momentum equation: 


de , . e du 1 

— + {2 + H)* - = -Cf 

dj] ^ U drj 2 ^ 


3.3 

The momentum integral equation is numerically integrated along a surface streamline. 

The laminar region of the boundary layer is modeled by numerically integrating 
the following exact differential equation. The equation is solved iteratively through 
numerical integration along a streamline starting at a stagnation point [Ref 5]; 



0.45v 

uinY 


|(l + 2.222g{K, fip{r,y dr, + 0(0)’ 

0 




3.4 


where; U - velocity at outer edge of boundary layer 

0 - momentum thickness 

V di] 

r) - generalized coordinate along a streamline 

The value g(K,u) is based on exact solutions for a number of pressure 
distributions. Initial work was conducted by Thwaites with improvements by Curie 
[Ref 10]: 

= F,{K) - juG,{K) - 0.45 6K 3.5 


CMARC uses an empirical transition model based on the average pressure 
gradient, K , for predicting laminar to turbulent transition. The following relations are 
used to calculate the transition point [Ref 5]: 
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3.6 


\Kdri 

K=— - 

n-n,ns 

where is the streamline coordinate at instability. And, K is the local pressure 
gradient at boundary layer instability [Ref. 5]: 

K = -0.4709 + 0.11066* ln(Re J + 0.0058591 * ln^(ReJ (o < Re^ < 650) 

3.7 

K = 0.69412 - 0.23992 *ln(Re^) + 0.0205 * In-(Re^) (650 < Re^ < lOOOO) 

The local Reynolds number at transition is correlated to K with the following expressions 
[Ref 5]: 


K = -0.0925 + 0.00007 * Re^ 

K = -0.12571 + 0.000114286 * Re^ 

K = 1.59381 - 0.45543 * In(Re^) + 0.032534 * In' (Re^) 


(O < Re^ < 750) 
(750<Re^ <1100) 3.8 

(ll00<Re^ <3000) 


At transition, the initial turbulent shape factor, H, is given by the following empirical 
formula that is a fit to data developed by Coles [Ref 10]; 


H = 


1.4754 

log,o(Re^) 


+ 0.9698 


3.9 


Provisions are made to check for turbulent reattachment if laminar separation is 
encountered. At laminar separation, a point calculation is made to determine if the 
boundary layer will reattach. If reattachment is predicted, the boundary layer code 
immediately switches to turbulent calculations. No attempt is made to model the laminar 
separation bubble or provide a transition length. After laminar separation is predicted, 
the following empirical relations are used to determine if reattachment occurs [Ref 5]; 
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(Re, >125) 
(Re, <125) 


3.10 


K = 0.0227 - 0.007575 * Re^ - 0.000001157 * Re/ 
K = -0.09 


The boundary layer code in CM 2 \RC uses a point transition model. No attempt is 
made to model a more representative transition length. Turbulent calculations begin at 
transition using the Nash-Hicks model [Ref. 5]. Calculations continue along the 
streamline until turbulent separation is predicted or the end of the streamline is reached. 
No boundary layer data is available after separation. 

The authors of PMARC caveat that their boundary layer calculations are quite 
accurate for predominately 2D flow but break down in regions of large cross flow near 
separation. This premise was tested in Reference [1] by comparing the predominately 2D 
flow over the inboard region of a high aspect ratio wing to the finite difference 
calculations performed by the Naval Postgraduate School Unsteady Potential Flow Code 
(UPOT). In general, CMi^RC provided correct trends for both the transition and 
separation points. However, in all cases, CMARC predicted early transition and late flow 
separation. As expected, the differences were greatest at lower Reynolds numbers where 
boundary layer thickness is larger. 

Reference [1] also modeled flow over an inclined prolate spheroid. The 6:1 
prolate spheroid was chosen because extensive experimental data is available. In 
addition, the three-dimensional flow around a prolate shperoid is similar to a streamlined 
slender fuselage. With proper wake placement, CMARC was found to produce accurate 
normal force and pitching moment coefficients. Over the three-dimensional body, 
CMiARC boundary layer calculations also predicted early transition and late flow 
separation. Despite inaccuracies, CMARC boundary layer calculations remained useful 
when used as a design tool for visualizing the trend in transition and separation points 
with configuration changes. 
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IV. AERODYNAMIC MODEL OF THE FROG UAV 


A. BACKGROUND 

The Naval Postgraduate School Aeronautics Department is integrating UAV 
hardware and software to demonstrate autonomous flight, trajectory tracking and 
automatic landing. A core requirement for flight control law development is a valid 
aerodynamic truth model for the UAV airframe. A panel code model of the FROG UAV 
is one method for estimating many of the stability derivatives required for an 
aerodynamic truth model. This development effort concentrates on finding the static, rate 
damping and control power derivatives for both the longitudinal and lateral-directional 
axes. 

Panel code modeling usefulness goes beyond the development of aerodynamic 
coefficients. Flight control systems require accurate pitot-static and angle-of-attack 
sensor inputs. CMARC accurately solves on-body static pressure distributions and off- 
body flow velocities over the predominately attached flow fields of fuselage fore bodies. 
In this study, correction curves are generated for static source and angle-of-attack probe 
position errors. 


B. FROG UAV DESCRIPTION 

The FROG UAV is a small single engine flight test vehicle used for autonomous 
flight research by the Naval Postgraduate School Aeronautics Department. The aircraft 
was originally designated the FOG-R by the U. S. Army. It was designed as a small 
lightweight, battlefield observation platform that could be guided by a fiber optic data 
link. Table 4.1 presents the basic aircraft specifications. 

The aircraft is somewhat unconventional in that the engine is mounted in a nacelle 
tractor style above the fuselage and wing. The aft fuselage consists of a 1.75 inch 
diameter aluminum tube which connects the tail surfaces to the main fuselage. Figure 4.1 
displays a three view drawing of the FROG UAV. 
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PARAMETER 

MEASUREMENT/UNITS 

Length 

8.125 ft 

97.5 in 

Height 

1.75 ft 

21 in 

Weight 

67.7 lbs 

Power Plant 

12 Hp/ 2 Cycle 

Wing Airfoil 

NACA 2415 

Horiz. Stab. Airfoil 

NACA 0006 (Approx.) 

Sw (Sref) 

17.57 ft^ 

2530 in^ 

St 

3.174 ft^ 

457.1 in^ 

Sv 

0.9818 ft^ 

141.4 in^ 

C 

1.66 ft 

20 in 

Ct 

0.958 ft 

11.5 in 

bw 

10.54 ft 

126.5 in 

bt 

3.313 ft 

39.75 in 

bv 

1.25 ft 

15.0 in 

It 

4.44 ft 

53.25 in 

Iv 

4.44 ft 

53.25 in 

ARw 

6.32 

ARt 

3.46 

ARv 

1.59 

Vh 

0.49 

Vv 

0.16 


Table 4.1 FROG UAV Characteristics, after Ref. [2]. 
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Figure 4.1 FROG UAV Three-View Drawing. 
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The FROG UAV, as operated by NFS, is equipped with airspeed, angle-of-attack, 
altitude and control surface sensors. In addition, a miniature Inertial Measurement Unit 
(IMU) captures aircraft attitude, acceleration and body rates. Data is down linked to a 
mobile SGI workstation through a spread spectrum modem. Onboard GPS provides 
differential GPS navigation capability with the ground station used as a reference. The 
aircraft can be flown by conventional radio control or by up-linking flight control 
commands from the SGI workstation. 

Current flight control development requires a linearized aerodynamic model 
around the cruise trim point of 60 m.p.h. or 88 ft/s. This flight condition is selected for 
the development of stability derivative data with CMARC. Table 4.2 lists the aircraft 
parameters for the trim flight condition. 


PARAMETER 

MEASUREMENT 

UNITS 

Weight 

67.73 

lbs 

IXX 

12.52 

slug-ft^ 

lYY 

8.43 

slug-ft^ 

IZZ 

18.55 

slug-ft^ 

Airspeed 

60/88 

mph and ft/s 

Altitude 

800 

ftMSL 

Air Density 

0.002327 

Slug/ft3 

Center of Gravity 

34.5% 

M.A.C 

Cl trim 

0.4295 

n/a 

^ trim (est) 

-0.04 

degrees 

^ Etrim 

5.1 

degrees 


Table 4.2 FROG UAV Trim Flight Condition, after Ref.[2]. 

C. FROG UAV MODELING 

1. General 

LOFTSMAN is utilized for the creation of all CMARC input file patches except 
for wing tips. In some cases, CMARC’s more efficient built-in capability to model 
standard NACA 4-digit wing surfaces could be used. However, with the requirement to 
mesh deflected control surfaces, all patches are created with LOFTSMAN from the start. 
The complete FROG UAV model with all patches activated is displayed in Figure 4.2. 
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Appendix A contains detailed descriptions of obtaining stability derivative data in a less 
formal user guide format. 



Figure 4.2 FROG UAV Panel Code Model with All Patches Activated. 

Some assumptions are made to simplify the modeling process. First, the 
horizontal and vertical stabilizers are modeled with a NACA 0006 section. The actual 
surfaces are constructed with a flat section, rounded at the nose and tapered starting at the 
control surface hinge line to a sharp trailing edge. The NACA 0006 provides a close 
approximation and allows the use of LOFTMAN’s built-in wing lofting capability. For a 
potential flow solution, this simplification is considered minor. 

A second simplification is made regarding the vertical stabilizer’s tip rib 
orientation. The actual rib is canted down 5® with respect to the longitudinal waterline. 
LOFTSMAN will only model a chord line that is parallel to the waterline (constant BL). 
Thus, for this study, the vertical tail tip rib is modeled with a constant BL, but the span is 
adjusted to maintain the same overall surface area. 

A third assumption is made regarding aileron control surface hinge lines. The 
actual FROG UAV ailerons are piano hinged along the upper surface. For simplification 
of the LOFTSMAN model, all surfaces are hinged at the centerline. 
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Finally, there is no attempt to model the tricycle landing gear struts or wheel 
assemblies. The landing gear components do not contribute significantly to the 
aerodynamic stability derivatives. However, they certainly need to be taken into account 
when measuring moments of inertia for a dynamic model. 

2. Modeling Coordinate System 

The model is developed using a coordinate system selected to simplify fuselage 
measurements. The +x-axis starts even with the nose and runs aft along the bottom of the 
fuselage, parallel with the tail boom. The bottom of the fuselage is used as the waterline 
with +z-axis in the up direction. This allows for easy vertical measurements when the 
aircraft is placed flat on a horizontal surface. The +y-axis runs fi'om centerline outboard 
parallel to the right wing. Figure 4.3, which displays static-pressure source and alpha 
vane locations, also shows the location and origin of the modeling coordinate system. 

3. LOFTSMAN Patches 

LOFTSMAN is used to generate all the model patches except for wing tips. 
CMARC’s built-in capability is used to create wing tip patches. Appendix B contains 
listings of all the LOFTSMAN input files. Once a surface is meshed, the mesh is saved 
to a file as a CMARC/PMARC patch. The resulting text file is then opened, and the text 
is copied and pasted with any text editor into the patch definition section of the CMARC 
input file. LOFTSMAN patch files are not listed because they are redundant with the 
patches in the final CMARC input file listed in Appendix C. Multiple patches for the 
same control surface may be stored in the input file. For instance, two horizontal 
stabilizer patches with different elevator deflection angles may be kept in the CMARC 
input file. However, only one may reside in the active patch input section. All other 
patches for the same aircraft component must reside after the patch with the "last patch" 
toggle set (TNODS=5). 

When saving a patch, LOFTSMAN automatically takes care of all CMARC input 
formatting. A patch, as formatted by LOFTSMAN, assumes additional patches will 
follow in the CMARC input file. Therefore, the last segment’s TNODS variable is set 
TNODS=3. When the patch is the last patch in the input file, the TNODS variable must 
be manually set to TNODS=5. If CMARC hangs up while reading in geometry 
information, most likely TNODS=5 is missing on the last patch. 


22 


a. Fuselage Model 

The fuselage is lofted as a B-type body. A B-type body is used when 
major portions of the fuselage have a circular or oval cross section. The input file is 
listed in Appendix B. Only the right side is meshed, with a symmetric left side created 
by toggling the IPATSYM variable to IPATSYM=1. LOFTSMAN assumes that B-type 
bodies converge to a specific point at the fore and aft ends. The flat aft fuselage face does 
not provide this single point. A slight modification was made to the aft face to allow 
automatic meshing as a B-type body. The center of the aft face is extended very slightly, 
approximately 1/8 inch, to provide a convergence point for the final rear triangular 
panels. This small deviation is assumed not affect the aerodynamic fidelity of the model 
for a potential flow solution. 

The right side was originally meshed separately fi-om the wing as a 20 x 
20 panel patch. This created a low order fit when the wing patch was butted to the side of 
the fuselage, resulting in overlapping panels. A final mesh was created that flowed 
around the wing root and fuselage intersection for a high order fit. All the fuselage 
panels at the wing root join with the adjacent wing panels. This mesh requires that the 
fuselage be broken up into six separate panels per side. They are the nose patch, the 
forward transition patch, the top and bottom wing root patches, the aft transition patch 
and finally the rear fuselage patch. Some manual editing is required to straighten out 
panels on the upper fuselage patch. When the six patches are added together, the final 
configuration is modeled with a 44x15 panel patch. 

b. Main Wing Patch 

The NACA 2415 wing is created with four separate patches to allow the 
addition of a deflected aileron mesh. CMARC comes with a broad selection of “*.SD” 
airfoil template files that are automatically loaded during installation. The 
“NACA2415.SD” file is used for this model. The 10x15 inboard wing patch runs fi’om 
the wing root, past the flaps, to the start of the aileron. The 10x15 mid-wing patch covers 
the portion of the wing spaimed by the aileron. The tapered outboard wing extension is 
made with a 6x15 patch. Finally, a 4x15 semi-circular wing tip patch is created using 
CMARC’s built-in wing tip functionality. The wing is set to a 4.5° incidence in the 
LOFTSMAN input file. Alternatively, the patch could be created with no incidence and 
then the patch coordinate system could be rotated in the CMARC input file. Together, 
the three wing patches add to make a 30x15 panel wing model. 
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The ailerons are meshed at zero degree deflection and five degrees of right 
rolling moment deflection. For the no deflection case, a single right wing patch is 
meshed with CMARC's symmetric patch toggle (IPATSYM=1) creating the left patch. 
For deflected ailerons, separate patches are meshed for each wing. The right wing aileron 
deflection is five degrees trailing edge up (T.E.U.) and the left wing five degrees trailing 
edge down (T.E.D.). 

c. Horizontal Stabilizer Patch 

The horizontal stabilizer and elevator patch is created with a single 10x16 
mesh using the “NACA0006.SD” airfoil template. As stated earlier, the actual UAV 
horizontal stabilizer has flat upper and lower surfaces and a roimded leading edge. The 
NACA 0006 airfoil is substituted because it provides a close approximation to the leading 
edged radius and thickness ratio. No other special modifications are required. A tip 
patch is not added because some of the resulting panels would be too small. In particular, 
the triangular panels closing out the aft end of the tip are too small in proportion to the 
other panels, creating singularities due to machine resolution. Initially, an attempt was 
made to model horizontal and vertical stabilizer wing tips, but the model would not 
converge with them. Leaving off tip patches does not significantly influence results 
according to the CMARC User's Guide [Ref. 3]. 

The elevator spans the entire horizontal stabilizer. Two patches are 
meshed, one with zero deflection and the second with positive five degrees (+5° T.E.D.) 
of deflection. The deflected patch is used to obtain the elevator control power derivatives 
(Cue> CiviSe and CoSe)- 

d. Vertical Stabilizer Patch 

Two different versions of the vertical stabilizer are used for FROG 
modeling. An extended vertical stabilizer, which includes the "effective" area of the 
empennage tail boom, is used to measure the static sideslip and yaw rate damping 
derivatives. It is a 10x14 panel patch. The second patch models the actual area of the 
vertical stabilizer. It is a 10x12 patch and is meshed with a five-degree rudder deflection 
for measuring control power derivatives. 

The vertical stabilizer and rudder patches are created with a single mesh 
using the “NACA0006.SD” airfoil template. As with the horizontal tail stabilizer, the 
NACA 0006 airfoil closely approximates the vertical surface leading edge radius and 
thickness ratio. The LOFTSMAN input file is different in that a vertical wing surface 
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requires a modification to the rib axis. The rib axis must be specified with an x-axis 
rotation of 90°, a y-axis rotation of 0° and an unspecified (999.0) z-axis rotation. No 
symmetry is selected for the vertical stabilizer because the patch is already S3nmnetric 
about the y=0 plane. As with the horizontal stabilizer, a tip patch is not added because 
some of the resulting panels would be too small. 

The rudder spans the entire vertical stabilizer. Two patches are meshed, 
one with zero deflection and the second with positive five degrees (+5° T.E.L.) of 
deflection. The deflected patch is used to obtain the elevator control power derivatives 
(Cu 6) CmSc Slid CoSe)' 

e. Tail Boom Patch 

The tail boom patch is created as a single 12 x 10 mesh using a B-type 
body. Again, only the right side is meshed due to symmetry. The LOFTSMAN input file 
requires modifications at both ends in a similar fashion to the aft fuselage. A single point 
is added to allow convergence of the triangular panels at either end. With this point, the 
tail boom has the appearance of being tapered at both ends. The point is then manually 
edited out in the CMARC input file by replacing the “x” coordinate of the beginning and 
ending section panels with the correct value. In most cases, the tail boom is left out of 
solution to aid in convergence. This is due to the small overlapping panels at the fuselage 
tail boom junction. Being a slender, round tube directly in the fuselage slipstream, the 
tail boom is assumed to have a small influence on the stability derivatives. 

f. Engine Pod Patch 

The engine pod patch, or nacelle, is created as a single 15 x 10 mesh using 
a B-type body. Only the right side is meshed due to symmetry. The prop spinner is an 
integral part of the patch. No attempt is made to model the prop, engine heads or exhaust 
system. As with the tail boom, in most cases the engine pod is left off the model used to 
establish stability derivatives. 

g. Engine Pylon Patch 

The engine pylon patch is modeled with a single 15x10 mesh using an A- 
type body. A-type bodies are used to model surfaces similar to boat hulls with cornered 
surfaces or sharp chines. In addition, A-type bodies do not require the body to be 
completely enclosed. As a result, an A-body was selected to model just the sides of the 
pylon. Only the right side is meshed due to symmetry. A low order fit is achieved with 
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the adjacent fuselage and engine pod panels. This results in questionable pressure 
distributions. As a result, the pylon patch was turned off for most configurations. A 
future attempt could be made to create a high order fit between the other patches or to 
model the pylon as a flat surface. This will probably require manual editing of the 
intersecting patches. Even with a high order fit, a wake cannot be added to the trailing 
edge of the pylon. It would impact the vertical tail siuface 

4. Common CMARC Input File Errors 

The patches created in LOFTSMAN are assembled into a single CMARC input 
file with any text editor. The default input file, which comes with CMARC, or any old 
input file may be modified. There are many errors that will cause CMARC to hang up 
without an error message. The two most common errors are forgetting to designate the 
last patch and incorrectly numbering the wake patches. 

The last patch must be designated by including a TNODS=5 setting in the last 
section of the last patch. If it is not included, CMARC hangs up when reading in the 
geometry. In a similar manner, the last wake must be designated with a NODEW=5 
setting. If the last wake is not designated, CMARC hangs up while reading in the wake 
information. Another common error involves incorrect wake to patch number 
association. Patch numbering changes whenever patches are disabled or reordered. The 
KWPACH field for each wake definition must be checked to make sure it reflects the 
current patch numbering. 

D. STATIC-PRESSURE SOURCE AND ALPHA VANE CORRECTIONS 

THROUGH OFF-BODY FLOW ANALYSIS 

CMARC is ideally suited for off-body flow analysis. Off-body streamlines may 
be placed through a point anywhere in the flow field. CMARC will then follow the 
streamline up and downstream the distance designated in the input file. This is 
particularly useful for flow visualization. In addition, CMARC calculates pressure 
coefficient and velocity at each point along the streamline. For this study, two 
streamlines are placed through the locations of the static-pressure source and alpha probe 
locations. Pressure coefficient is used to quantify static source position error and velocity 
is used to calculate alpha probe position error as a function of FROG UAV angle-of- 
attack. Both static pressure and AOA are digitized for down link to the ground station 
allowing the values to be easily corrected. Either a look-up table or curve fit correction 
can be applied subsequent to being passed to the flight control routines. This analysis 
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was previously published by the author in Ref. [1], but is included here to provide a 
single source of UAV modeling conducted to date. 

1. Description of the FROG UAV Pitot-Static and AOA Systems 

The pitot-static system and angle-of-attack probe share a common flight test 
boom extending from the nose of the UAV. The boom contains both the total and static 
pressure ports. Figure 4.3 depicts the general dimensions of the flight test boom 
installation and the modeling coordinate system. 

2. Modeling Off-Body Streamlines 

Streamlines are placed at the two locations indicated in Figure 4.3, which 
correspond to the static source and alpha probe locations. Two off body streamlines were 
activated in CMARC by setting NSTLIN=2 in the &SLIN1 line. Only a short distance of 
2 inches is selected up and downstream in the SU and SD fields to reduce the size of the 
output file. Figure 4.4 is a POSTMARC rendering of the two off-body streamlines used 
for sensor corrections. With the model at at=0°, notice that the streamline is curving up 
at the angle-of-attack vane location 6.5 inches in front of the aircraft nose. 

3. Analysis of Static Source Position Errors 

The position error pressure coefficient, AC? pc or APp/qc, is a fimction of free 
stream Mach number and angle-of-attack provided that the static source is located outside 
of a thick boundary layer and sideslip is minimized [Ref 11]. In the case of the FROG 
UAV with incompressible flow, APp/qc becomes a fimction of angle-of-attack only. 
Therefore, corrections can be simply defined as a fimction of measured angle-of-attack. 

A DOS batch file was executed to step the CMARC model through angles-of- 
attack ranging from -8° to 20°. The batch file incremented the angle-of-attack using 
CMARC’s command line override feature. In addition, a new output file name was 
designated for each angle-of-attack. Position error pressure coefficient is then read from 
the off-body streamline listing of the output file at the location corresponding to the static 
source. Table 4.3 lists the values of APp/qc calculated from CMARC data. Figure 4.5 
displays ACp pc as a fimction of indicated angle-of-attack. The second order influence of 
angle-of-attack is clear with a second order curve fitting tightly through the data points. 
Of note, the error is relatively constant for a ±8° band around trim angle-of-attack. For 
incompressible flow, position error pressure coefficient is independent of airspeed and 
altitude. 
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Figure 4.4 FROG Off-Body Streamline Visualization with POSTMARC (a=10°). 

Position error pressure coefficient can be turned into position corrections for 
airspeed and altitude. The following relations were developed which assume small errors 
and incompressible flow; 

ViACp 

AVpc=—^ and AVp,=V,-Vi 4.1 

AVpcVj 

AHpe=—^ and AHpc=Hc-Hi 4.2 

^stdSo 

Where: 


AHpc is the altitude position correction. 

AVpc is the velocity position correction. 

APp 

ACn = —- or position error pressure coefficient, 
qc 

Cstd is standard day density ratio, 
go is the gravitational constant. 
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Table 4.3 displays corrections calculated for both airspeed and altitude at the 
FROG UAV trim condition of 88 fl/s and 800 ft MSL. The corrections are added to the 
indicated value to obtain the corrected value. Figures 4.6 and 4.7 display the corrections 
as a function of indicated angle-of-attack. Again, a second order curve fits nicely through 
the data points. Equations 4.1 and 4.2 can be used to implement a correction algorithm 
based on airspeed and altitude. 

4 . Analysis of Alpha Vane Position Error 

Local flow field velocity is extracted from the off-body streamline listing to 
obtain local angle-of-attack. The alpha vane is assumed to capture the x-z component of 
the local velocity field and ignore cross flow in the y direction. Flow field velocity is 
turned into indicated angle-of-attack and angle-of-attack position correction with the 
following equations: 





180 

71 


degrees 


4.3 


o 

Attpc =at -aj degrees 


4.4 


A DOS batch file is executed to step the CMARC model, with an off-body 
streamline located at the vane position, through angles-of-attack ranging from -8° to 20°. 
Local velocity components are then read from the location corresponding to the alpha 
vane. Table 4.4 lists the values of Attpc calculated from CMARC data. Figure 4.8 
displays Aa pc as a function of indicated angle-of-attack. Linear and second order curve 
fit equations are also indicated on Figure 4.8. Angle-of-attack correction is fairly linear 
through the FROG operating envelope, with approximately -1.25 degrees of position 
error at the FROG cruise trim condition. The corrections apply at all incompressible 
airspeeds and all altitudes. 

5. Summary of Off-Body Flow Field Analysis 

CMARC proved useful for both static-pressure source and alpha vane position 
corrections. Measured data may be corrected using look-up tables with the values in 
Table 4.3 and 4.4 or by using the curve fits in Figures 4.5 through 4.8. Flight testing is 
recommended for validation of sensor corrections obtained from this CMARC off-body 
flow field analysis. 
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UAV AOA 

aj (deg) 

ACppc 

AP/qc 

^ Correction 
AVpc=Vc-Vi(ft/s) 

H Correction 
AHpc=Hc-Hi (ft) 

-8 

0.1092 

4.8 

13.5 

-6 

0.1120 

4.9 

13.8 

-3 

0.1141 

5.0 

14.1 

-2 

0.1140 

5.0 

14.1 

-1 

0.1137 

5.0 

14.1 

0 

0.1132 

5.0 

14.0 

1 

0.1123 

5.0 

13.9 

2 

0.1111 

4.9 

13.7 

3 

0.1096 

4.8 

13.5 

4 

0.1078 

4.8 

13.3 

5 

0.1057 

4.7 

13.1 

6 

0.1034 

4.6 

12.8 

8 

0.0977 

4.3 

12.1 

10 

0.0909 

4.0 

11.2 

12 

0.0831 

3.7 

10.3 

14 

0.0741 

3.3 

9.2 

16 

0.0641 

2.8 

7.9 

18 

0.0530 

2.3 

6.6 

20 

0.0410 

1.8 

5.1 


Table 4.3 Position Error Corrections for the NPS FROG UAV at V=88 ft/s and 
H=800 ft MSL. Derived from CMARC Panel Code Off-Body Flow 
Field Analysis. 
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ACp pc 


Source; CMARC Panel Code 
Trim Condition; V=88 ft/s H=800 ft MSL 



AOA Indicated (deg) 


Figure 4.5 Position Error Pressure Coefficient, ACp pc for the NPS FROG UAV. 

Derived from CMARC Panel Code Off-Body Flovk^ Field Analysis. 
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AHpc (ft) 


Source. CMARC Panel Code 
Trim Condition V=88 tt/s H=800 ft MSL 
Flight Test Boom Static Source (8.25" forward of nose) 



AOA Indicated (deg) 


Figure 4.6 Altitude Position Error, AHpc, for the NPS FROG UAV at V=88 ft/s 
and H=800 ft MSL. Derived from CMARC Panel Code Off-Body 
Flow Field Analysis. 
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Source; CMARC Panel Code 
Trim Condition; V=88 ft/s H=800 ft MSL 
Flight Test Boom Static Source (8.25" forward of nose) 



AOA Indicated (deg) 


Figure 4.7 Airspeed Position Error, AVpc, for the NPS FROG UAV at V=88 ft/s 
and H=800 ft MSL. Derived from CMARC Panel Code Off-Body 
Flow Field Analysis. 
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UAV AOA 

at (deg) 

Velocity at Alpha Vane 

AOA Correction 
Aa=aT-<x, (deg) 

AOA Indicated 

tti (deg) 

Vx (ft/s) 

Vy (ft/s) 

Vz (ft/s) 

-8 

80.92 

1.66 

-12.23 

0.60 

-8.60 

-6 

81.27 

1.65 

-8.65 

0.08 

-6.08 

-3 

81.60 

1.64 

-3.21 

-0.75 

-2.25 

-2 

81.67 

1.63 

-1.47 

-0.97 

-1.03 

-1 

81.71 

1.63 

0.28 

-1.20 

0.20 

0 

81.73 

1.62 

2.13 

-1.49 

1.49 

1 

81.73 

1.61 

3.93 

-1.75 

2.75 

2 

81.70 

1.60 

5.72 

-2.00 

4.00 

3 

81.66 

1.59 

7.51 

-2.25 

5.25 

4 

81.58 

1.58 

9.30 

-2.50 

6.50 

5 

81.48 

1.57 

11.08 

-2.75 

7.75 

6 

81.37 

1.56 

12.88 

-2.99 

8.99 

8 

81.07 

1.53 

16.43 

-3.46 

11.46 

10 

80.67 

1.51 

19.98 

-3.91 

13.91 

12 

80.17 

1.48 

23.50 

-4.34 

16.34 

14 

79.61 

1.46 

26.99 

-4.73 

18.73 

16 

78.93 

1.43 

30.47 

-5.11 

21.11 

18 

78.18 

1.39 

33.90 

-5.44 

23.44 

20 

77.34 

1.36 

37.31 

-5.75 

25.75 


Table 4.4 Angle-of Attack Vane Position Error Corrections for the NPS FROG 
UAV. Derived from CMARC Panel Code Off-Body Flow Field 
Analysis. 
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Source: CMARC Panel Code 
Trim Condition: \/=8d ft/s H=800 ft MSL 
Flight Test Boom Alpha Vane (6.5" forward ol nose) 



Figure 4.8 Angle-of-Attack Vane Position Error, Attpc, for the NPS FROG UAV. 

Derived from CMARC Panel Code Off-Body Flow Field Analysis. 
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E. 


DEVELOPMENT OF STABILITY DERIVATIVES 


In this section, CMARC is used to develop the longitudinal and lateral-directional 
stability derivatives for the FROG UAV. The development effort focuses on finding the 
static, rate-damping and control-power derivatives. The results obtained from CMARC 
are compared to data sets obtained from empirical estimation techniques and flight-test 
parameter estimation. In the next section, the stability derivatives are entered into a 
dynamic model to find modal frequency, damping and response to control deflections. 

The potential flow analysis performed by CMARC does not provide accurate 
viscous drag values. Therefore, total drag is estimated from the flight-test power-off 
glide ratio and cruise thrust required. The User Guide developed in Appendix A 
describes in detail how to gather stability derivative data with CMARC. An abbreviated 
discussion is presented below. 

CMARC contains built-in frmctionality to integrate forces and moments in all 
axes over the surface of a body. Force and moment coefficients are automatically 
nondimensionalized based on the mean aerodynamic chord, reference wing area, semi¬ 
span and center of gravity location in the CMARC BINP9 input line. Coefficients are 
presented in both wind and body axes. Of note, CMARC uses the semi-span to 
nondimensionalize rolling and yawing moment coefficients. Most texts on stability and 
control, including Roskam [Ref. 12] and Etkin [Ref 13], nondimensionalize rolling and 
yawing moments by span. Rolling and yawing moments are nondimensionalized by span 
in this study. Table 4.5 summarizes the factors for normalizing moments and angular 
rates. All rolling and yawing moment coefficients presented in this study have been 
normalized with span by dividing the CMARC output by a factor of two. Table 4.5 also 
indicates the characteristic time, t*, employed for angle rate data reduction. 

In addition to differences in nondimensionalizing moments, CMARC uses the 
typical CFD axis system shown in Figure 4.9. For this study, all work is performed in the 
stability axis system. Figure 4.9 also illustrates the standard stability axes implemented 
in this study. The sign of CMARC roll and yaw moments need to be reversed. The 
direction for positive control deflections is also shovm in Figure 4.9. All control surfaces 
are patched with positive deflections using the convention in Figure 4.9. 
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A potential flow solution will not produce satisfactory results for bodies with 
significant areas of flow separation. Therefore, CMARC models must be analyzed in the 
linear slope regions for valid results. 


MOMENTS 

NORMALIZING 

PARAMETER’ 

RATES 

CHARACTERISTIC 

TIME 

L = CiqSb 

b 

II 


^ = CmQSc 

c 

rc 

r = - 

2u^ 


N = C,qSb 

b 

rb 

r- - 

» 

II 

S 

o 


Note: 1) CMARC nondimensionalizes roll and yaw coefficients with b/2. 


Table 4.5 Nondimensionalized Moment and Rate Equations. 

For the static stability derivatives, the CMARC model is run at two different 
angles-of-attack and one sideslip angle with zero control surface deflections. The slopes 
of the force and moment coefficients are then taken to produce the Cta and Cma 
longitudinal derivatives and the Cyp, Cip and Cnp lateral-directional derivatives. 

For the control-power derivatives, the model is run at the trim condition with 
successive control deflections. The difference between the results with and without 
control surface deflections yield the Ctse, Cmsb and Cose longitudinal and the Cysr, Cisr, 
CnSr, Ci5a and CnSa lateral-directional control-power derivatives. 

Development of the damping derivatives is not as straight forward. The static 
derivatives are obtained with motion disabled. For the longitudinal damping derivatives, 
the model is run with oscillating vertical plunging motion to obtain the Cl and Cm a-dot 
terms. The lift and pitching moment coefficients are broken into real (in phase with 
AOA) and imaginary (out of phase with AO A) components. The imaginary components 
are due to a-dot effects. Next, the model is run with oscillating pitch motion to obtain 
the combined a-dot and pitch rate terms. Subtracting the a-dot influence obtained from 
the plunging motion isolates the pitch rate-damping term from the pitch motion. 
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CMARC Axis System 



Stability Axis System 


Figure 4.9 CMARC CFD Axis System Compared To Stability Axis System. From 
References [3] and [14] Respectively. 
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For the lateral-directional analysis, the P-dot terms are assumed to be small. This 
allows the model to be run with just oscillating roll and yaw motion. As with the 
longitudinal test case, the imaginary or out of phase component yields the combined P- 
dot and rate-damping data. With the p-dot terms assumed to be small, the oscillating 
motion yields the Cy, Ci, and Cn roll and yaw rate terms directly. 

Throughout this chapter, stability derivatives obtained from CMARC are 
compared to three other data sets. The first comes from classical analysis through the 
methods presented by Roskam in References [12] and [15]. The second set comes from a 
classical aerodynamic model developed by Papageorgio [Ref 2] and further refined 
through parameter estimation by Engdahl. The resulting parameter estimation model 
provides a good match to observed flight characteristics and is currently in use by the 
NPS Department for Aeronautical Engineering for closed-loop flight control 
development. The third set of data is for the Cessna 172 at the cruise flight condition. 
The data is only presented for an order of magnitude comparison. 

1. Longitudinal Stability Derivatives 


Static Longitudinal Stability Derivatives 


Three basic longitudinal stability derivatives can be measured with just 


two runs of the CMARC model. The model is first analyzed at an angle-of-attack 
corresponding to the estimated trim condition. In this case, at=0° is selected for the first 
run. A second CMARC run is conducted with angle-of-attack incremented by one or two 
degrees. Cl and Cm are then extracted manually from the data files. The slope of Cl and 
Cm versus angle-of-attack provide the Clo and Cma longitudinal derivatives. For this 
study, several angles-of-attack were analyzed to check consistency of the slope. In 
addition, amm is calculated from the lift curve slope and trim lift coefficient. For the 
longitudinal analysis, only half the model is analyzed. The symmetric calculation mode 
is selected by setting both RSYM=0.0 and IPATSYM=0 in the CMARC input file. 
Equations 4.5 through 4.7 are used for these calculations: 


^ (Cl,-Cl,) . 180 

Lrf 



per radian 


4.5 
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per radian 


4.6 


C 


ma 


(Cm2 Cmi) ^ 18Q 

(a2-ai) 71 


a 


trim 


= a°, 



180 

71 


degrees 


4.7 


Two FROG UAV model configurations are analyzed in a build-up 
approach to check results against classical calculations and flight-test data. Figure 4.10 
shows the CMARC models. First, just the wing and horizontal tail are considered. The 
patches for all other surfaces and wakes are turned off and the wing root is extended to 
centerline. Then, the blended wing/fuselage and horizontal tail are analyzed. Values of 
Clo and Cma for these two configurations are presented in Table 4.6. 

Classical design calculations are also performed to estimate Cma for 
comparison to CMARC results. Equation 4.8 is used for the calculation of Cma^ 




(h-h,,)-VH-^ 1- 


w 


ds 

daJ 


4.8 


In classical design, the horizontal tail downwash derivative, ds/da, is generally selected 
fi'om empirical data. Using a taper ratio of TR=1:1 and an aspect ratio of AR=6, 
ds/da=0.4 is selected from empirical charts in Ref [16] for the FROG UAV 
configuration. Also, calculations for ds/da=0.25 are included to illustrate pitching 
moment sensitivity to the downwash derivative. Table 4.6 lists Cma for oach 
configuration. 

The final static longitudinal parameter required is total aircraft drag. Drag 
coefficient plays an important role in long-period aircraft dynamics. Unfortunately, 
potential flow panel codes such as CMARC do not provide accurate total drag estimates. 
They can provide good induced drag predictions. And, if equipped with a boundary layer 
code like that contained in CMARC, they can provide integrated skin fnction results. 
However, a large total drag contribution in the form of pressure drag due to flow 
separation is not accounted for. Total drag estimates are made below using the two 
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simple techniques shown in Equations 4.9 to 4.12. The first method is based on flight- 
test glide ratio. The second is based on cruise power required and estimated prop 
efficiency. Note that the selected prop efficiency is relatively low due to the small 
propeller diameter, high RPM and pusher configuration. The two methods provide drag 
predictions within 10% of each other. The results are averaged to Cd=0.065 and included 
in Table 4.6. 

Method 1; Lift-to-Drag Ratio (L/D=7 from flight-test) 



4.9 


D 


9.67 lbs 


= 0.0611 


0.5*0.002327 lb ■ ft^ls^*\1.51 ft 


4.10 


Method 2; Cruise Power Setting (HP=5, rip=0.35) 



D HP„*550 ftlbslslHP*'nplV 

qS qS qS 


qS 


4.11 


5 HP* 550 ft-lbs/s/HP*0.35/SS, ft/s 


4.12 


0.5*0.002327 lb-s^/ff*SS^ /rV5^ *17.57 ft^ 


Average: Cpave =(0.0611+0.069)/2 = 0.065 
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Simplified Wing and Horizontal Tail 



Blended Wing and Fuselage with Horizontal Tail 


Figure 4.10 CMARC Models of the FROG UAV for the Determination of Static 
Longitudinal Stability Derivatives. 
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Flight-test data for the short-period and phugoid modes were used for 
longitudinal parameter estimation. Values for Cta ^d Cma based on preliminary 
parameter estimation work by Engdahl are presented in Table 4.6. 


METHOD 

CONFIGURATION’ 

STATIC 

^rim 

{oeQ) 

: LONGIT 

Clo 

^per rao) 

UDINAL 

^ma 
(per rao) 

PARAME 

Cq 

TERS 

Cdo 
iper rao) 

CMARC 

Wing/Horiz Tail 

-0.87 

4.86 

-0.835 

n/a 

n/a 

Panel Code 

Wing/Fuselage/Horiz Tail 

-0.01 

4.85 

-0.413 

n/a 

0.266 

Classical 

Wing/Horiz Tail - 5e/5a=0.25 

-0.81 

4.85 

-1.00 

n/a 

n/a 

Design^ 

Wing/Horiz Tail - 58/5a=0.40 

-0.82 

4.82 

-0.70 

n/a 

0.253 

Parameter 

Estimation'* 

Flying Aircraft 

n/a 

4.09 

-0.42 

0.0655 

n/a 


C-172® 

n/a 

0.31 

-0.89 

0.31 

0.13 


NOTES: 1) CGx=34.5% M.A.C. / CG2=8.6” from bottom of fuselage. 

2) Zero lift wing incidence is +6.5° from the longitudinal reference line. 

3) Classical design after Ref. [12]. 

4) Parameter estimation from NPS flight test data by Engdahl. 

5) Average drag obtained from L/D and cruise power analysis from flight test data. 

6) C-172 data from Ref. [12]. 

Table 4.6 Comparison of FROG UAV Static Longitudinal Stability Derivatives. 

Good correlation of lift curve slope (Cta) to classical design is 
demonstrated by the two CMARC configurations in Table 4.6. Both techniques over 
estimated Clo compared to flight-test data. Neither the classical design technique nor the 
CMARC model accounts for control surface leakage or separation effects. Modeling 
actual flight control gap geometry may provide a closer match to experimental data. 

Excellent correlation is achieved for Cma between the blended 
wing/fuselage CMARC model and flight-test results. The low value of Cma implies that 
ds/da is greater than the 0.4 value selected from empirical charts in Reference [16]. The 
simplified wing/tail model resulted in a Cma that closely matches classical design 
predictions with d£/da=0.35. Clearly, CMARC modeling provides the value-added 
capability to capture fuselage influences that are missed with classical design methods. 
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An investigation was performed to see if improved wake definition could 
better capture the ds/da downwash derivative. Initially a flexible wake was selected. 
However, at the cruise angle-of-attack, the flexible wake impacted the horizontal tail 
causing inaccurate results. A further investigation was made into a streamline based 
wake definition. A study by Walden et al. [Ref 17] studied wake turbulence by 
modeling an aircraft flying in trail of a wake-generating wing. A horizontal tail trailing 
the main wing is a similar configuration. The study found that a streamline-based wake 
is the best method for modeling downwash effects. In the study, the model was initially 
run with a rigid-wake. A streamline was also defined downstream of the main wing 
trailing edge. After analyzing the first results, the wake was then predefined to follow the 
streamline. The result was a first iteration on modeling a flexible wake with a predefined 
rigid-wake. This technique was tried during the course of FROG UAV modeling. 
However, FROG geometry proved unsuited to this technique. The streamline trailing the 
wing ended up impacting the horizontal tail. This indicated that a predefined wake would 
also impact the horizontal tail, resulting in inaccurate result. Figure 4.11 shows the 
streamline impacting the horizontal tail at cruise angle-of-attack with rigid-wakes 
selected. With the already close correlation of the rigid-wake results to flight-test data, 
further efforts to investigate Cma were abandoned. 



Figure 4.11 Rigid-Wake Streamline for Investigation of Streamline Based Wake 


Definition. 
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In summary, for the static longitudinal stability derivatives, CMARC 
produced accurate values for Cma a^id a slightly high value of Cua- Drag coefficient is 
obtained by averaging the values from the lift-to-drag ratio and cruise thrust required 
techniques. 


b. Longitudinal Damping Stability Derivatives 

The two aircraft motions illustrated in Figure 4.12 are used to develop the 
longitudinal rate-damping derivatives. A sinusoidal vertical plunging motion isolates the 
a-dot effects. Oscillatory motion is controlled with the CMARC BINPSA input file line. 
The oscillating pitch angle captures both the a-dot and pitch rate-damping terms. It is 
controlled with the CMARC BINP8B input file line. The a-dot influence developed 
from the plunging motion is subtracted from the pitching motion to isolate the pitch rate 
effect. All motion is conducted at a frequency of In rad/s, which equates to a reduced 
frequency of k=0.0595 for this configuration and trim airspeed. 



Figure 4.12 Aircraft Motion for the Determination of a-dot and Pitch Rate 
Damping. 
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The two models shown in Figure 4.13 are used for the rate-damping 
investigation. The complete wing/fuselage horizontal tail model is used for both the a- 
dot and pitch rate-damping terms, while the simplified horizontal tail only model is used 
to obtain a first approximation of pitch rate-damping. The User Guide in Appendix A 
provides a more detailed description of the techniques used to gather the dynamic 
longitudinal stability derivatives. 



Simpiifled Horizontal Tail Model 



z X 

Blended Wing and Fuselage with Horizontal Tail 

Figure 4.13 CMARC Models of the FROG UAV for the Determination of 
Longitudinal Rate-Damping. 
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Traditional design calculations, from Roskam [Ref. 12], are used for 
comparison to CMARC and flight-test results; 


f =2*C *n *V 


ds 

da 


4.13 


/ ^ 

C =-2*C *77 *V — 

^ ^ c da 


4.14 


p =9*C *T7 *V 

C =-1*0 *n *V — 

- ^Loh Ih ^ ^ 

where the dynamic pressure ratio at the tail is taken as Ph=1 0 

The dynamic longitudinal stability derivatives are presented in Table 4.7. 
CMARC-obtained values are compared to classical design values and C-172 stability 
derivatives. Flight-test data is unavailable for the damping derivatives. The aerodynamic 
model currently in use by the FROG research team implements longitudinal damping 
coefficients obtained from classical design techniques. The Cta-dot and Cmq obtained 
from CMARC for the complete model provided a good match to the classical design data. 
However, both Cta-dot and Cmq are over estimated by approximately 50% as compared 
to classical design values. One possible explanation for the additional damping from 
CMARC is that the classical approximation only estimates the horizontal tail influence on 
damping. CMARC also measures wing and fuselage damping influence. The horizontal 
tail by itself produces values close to the classical calculations. This supports the 
assertion that CMARC also captures the wing and fuselage influence for the whole 
model. 


4.15 


4.16 


48 


METHOD 

CONFIGURATION’ 

DYNAMIC LONGITUDINAL PARAMETERS 

C La-dot 
toeg; 

^Lq 

(per rao) 

'^ma-dot 

(per raa) 

^mq 
(per raoj 

Panel Code 

Blended Wing-Fuselage/Horiz Tail 

1.42 

6.82 

-6.24 

-11.78 

Horizontal Tail only 

n/a 

4.37 

n/a 

-11.94 

Classical 

Wing / Horizontal Tail 

1.56 

3.89 

-4.14 

-11.39 

Par. EsP 

Flying Aircraft 

n/a 

n/a 

n/a 

n/a 


C-172^ 

1.7 

3.9 

-5.2 

-12.4 


NOTES: 1) CGx=34.5% M.A.C. / CG 2 = 8 . 6 ” from bottom of fuselage. 

2) Classical design after Ref. [12], 

3) Parameter estimation from flight test data by Engdahl. 

4) C-172 data from Ref. [12]. 

Table 4.7 Comparison of FROG UAV Dynamic Longitudinal Stability Derivatives, 


c. Longitudinal Control-Power Derivatives 


The elevator control-power derivatives are obtained by substituting a 0° 
deflection horizontal-tail patch for one with +5° trailing edge down deflection. Only one 
run is required. The difference between the trim condition and the deflected value is 
divided by the elevator deflection as shown below: 


C, 




C, 


Dse 


-=-^*57.3 per rad 

^^^*2 o ^ 

- — -^*57.3 per rad 

Se^ 


4.17 


4.18 


C 




^Se 


5e2 -(^i 


■ 57.3 per rad 


4.19 


Elevator control-power derivatives are presented in Table 4.8. A classical 
design estimate of elevator control-power is also shown for comparison. The following 
relationships from Roskam [Ref. 15] are used for estimating elevator control-power: 
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4.20 


C, =C. 

ISp ^ 


C =-C *— or C =-C *V 


4.21 


where Clsf is the variation of lift coefficient with elevator deflection found from charts in 
Roskam [Ref 15]. 

In general, the elevator control-power derivatives obtained with CMARC 
correlate well with classical design techniques. For both Cue and CmSe, CMARC derived 
elevator power is approximately 15 percent higher than obtained from classical 
calculations. However, CMARC significantly under-estimates elevator control-power 
compared to the parameter estimation model. The most likely explanation for this 
phenomenon is that the CMARC solution does not model the effects of prop wash. The 
higher dynamic pressure over the horizontal stabilizer from prop wash increases elevator 
effectiveness. The parameter estimation model captures this phenomenon. This same 
trend can be seen in the modeling of rudder control-power. An attempt could be made to 
model the prop disk effects using CMARC's normal panel velocity definition capability. 


METHOD 

CONFIGURATION^ 

ELEVATOR CONTRO 

POWER 

Cuae 
(per raaj 

Cpse 
(per rao) 

^m6e 
(per raoj 

Panel Code 

Blended Wing-Fuselage/Horiz Tail 

0.438 

0.01 

-1.199 

Classical^ 

Wing / Horizontal Tail 

0.39 

n/a 

-1.04 

Par. Est^ 

Flying Aircraft 

1.13 

n/a 

-1.62 


C-172^ 

0.43 

0.06 

-1.28 


NOTES; 1) CGx=34.5% M.A.C. / CGz=8.6" from bottom of fuselage. 

2 ) Classical design after Ref. [12]. 

3) Parameter estimation from NPS flight test data by Engdahl. 

4) C-172 data from Ref. [12].. 

Table 4.8 Comparison of FROG UAV Elevator Control-Power Derivatives. 
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2. Lateral Directional Stability Derivatives 


a. Static Lateral-Directional Stability Derivatives 

Development of the static lateral-directional stability derivatives is not as 
straightforward. Both sides of the airframe must be modeled by setting both RSYM=1.0 
and 1PATSYM=1. This creates symmetric patches around the y=0 plane allowing 
CMARC to perform asymmetric calculations around the entire body and significantly 
increases processing times. 

Initially, a lateral-directional solution was attempted with a complete 
airframe model. The longitudinal model was modified by adding the vertical stabilizer. 
This required what was perceived to be a minor wake modification. Figure 4.14 
illustrates the configuration. The wing wake was terminated at the wing root to prevent 
interference with the vertical stabilizer. This left a wake gap for the vertical stabilizer. 
The engine nacelle and pylon patches were turned off because the wakes would also 
impact the vertical stabilizer. This configuration met with mixed success. The values of 
both Cyb and Cnb were opposite the expected directions. Under closer inspection, it was 
determined that the false vortex shed from the wing root caused a destabilizing influence 
on the vertical stabilizer. Figure 4.15 illustrates the vortex using CMARC's built in off- 
body streamline capability. 

In order to obtain satisfactory results, the model is broken up into two 
separate groups. The aircraft is modeled with the blended wing and fuselage as one 
group and the horizontal and vertical stabilizers as another group. Separate solutions are 
obtained for each group and summed through superposition to obtain the whole aircraft 
solution. The two separate models are shown in Figure 4.16. This method should be 
used as a last resort. A complete solution should be used if the airframe geometry allows 
adequate wake separation from the vertical stabilizer. A complete solution would capture 
fuselage, wing and vertical stabilizer interactions. 

Once the model is defined, it is checked for lateral-directional balance at 
zero sideslip angle (yaw angle=0°). The side force, rolling and yawing coefficients 
should be zero when a trial run is performed at zero sideslip. If lateral-directional forces 
or moments are present, the model and wake geometry should be checked for symmetry. 
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Figure 4.14 Unsuccessful Lateral-Directional Model of the FROG UAV. False 
Wing Root Vortex Caused Destabilizing Influence on the Vertical Stabilizer. 



Figure 4.15 Destabilizing Wing Root Vortex Visualized with Off-body Streamlines. 
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Figure 4.16 Final Simplified LateraFDirectional Models of the FROG UAV. 
Separate Solutions are Summed to Obtain Complete Airframe Solution. 

Next, a single CMARC run is performed with a= aoim and with two 
degrees of yaw angle. The results from both models are summed. The lateral-directional 
derivatives, Cyp, Cip and Cnp, are then obtained directly with equations 4.22 through 4.24: 


Cy = per radian 

c, 

Ap° 

r —Sjl* 

~ Ap° 


TV 

180 

n 

180 

71 


per radian 
per radian 


4.22 


4.23 


4.24 


It should be noted that the CMARC wind axis is modeled with x-aft and z- 
up, vice x-forward and z-down for the flight dynamics stability axis. Care must be taken 
to reverse the signs of the appropriate coefficients to convert from the CMARC wind axis 
to the flight dynamics stability axis system. 
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Static lateral-directional stability derivatives obtained from CMARC are 
presented in Table 4.9. For comparison three other sets of data are also presented. The 
first comes from classical analysis using methods from Roskam [Ref 12]. The second 
set comes from estimates based on data recorded from flight-test sideslip maneuvers, 
published by Papageorgio in Ref [2]. The third set comes from parameter estimation, by 
Engdahl, based on dynamic flight-test maneuvers. 


METHOD 

CONFIGURATION^ 

STATIC LAT-DIR PARAMETERS 

Cyp 

(per raa) 

(per raa) 

Cn(5 

(per raa) 

CMARC 

Panel Code 

Wing/Fuselage + HorzA/ert Stabs 

-0.249 

-0.063 

0.063 

Classical 

Design^ 

Wing/FuselageA/ert Tail 

-0.511 

-0.055 

0.051 

Flight Test^ 

Flying Aircraft 

-0.700 

-0.053 

0.057 

Parameter 

Estimation‘s 

Flying Aircraft 

-0.987 

-0.094 

0.176 


C-172® 

-0.310 

-0.089 

0.065 


NOTES; 1) CGx=34.5% M.A.C. / CG 2 = 8 . 6 " from bottom of fuselage. 

2) Classical design calculations by Roskam's methods, after Ref. [12], 

3) Flight test results from steady heading sideslip, from Ref. [2] 

4) Parameter estimation from flight test data by Engdahl. 

5) C-172 data from Ref. [12]. 

Table 4.9 Comparison of FROG UAV Static Lateral-Directional Stability 
Derivatives. 

CMARC produced weak results for the Cyp side force derivative. A large 
component of Cyp, approximately 60%, comes from the fuselage. Side forces are not 
modeled well by the potential flow solution. In addition, the engine pylon and pod are 
left off the model due to problems with their wakes impacting the vertical tail. They 
most likely provide a significant contribution to side force. 
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CMARC results for weathercock stability, Cnp, and dihedral effect, Cip, 
show close correlation to the classical calculations. However, dynamic parameter 
analysis indicates a considerably higher value for both derivatives. The parameter 
estimation most likely captures the additional dynamic pressure over the vertical 
stabilizer from prop wash. As with the longitudinal solutions, an attempt should be made 
to model the propeller disk with a normal velocity distribution. 

b. Lateral-Directional Damping Derivatives 

As with the static case, the dynamic solutions are run with separate 
wing/fuselage and vertical/horizontal stabilizer models. The solutions are summed 
through superposition. Figure 4.17 shows the two angular motions selected to develop 
the lateral-directional rate-damping derivatives. The (3-dot terms are generally 
considered to be small. As a result, the sideways plimging motion test case is not run. 
Yaw rate-damping is obtained directly from the oscillating yawing motion without the 
requirement to subtract (3-dot effects. Yawing motion is conducted at a frequency of Iti 
rad/s, which equates to a reduced frequency of k=0.369 for this configuration and trim 
airspeed. 





Roll Damping Motion 



Yaw Damping Motion 


Figure 4.17 Aircraft Motion for the Determination of Roll and Yaw damping. 
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The roll axis is de-coupled from both angle-of-attack and sideslip. 
Therefore, roll damping can be obtained directly from a constant rolling motion, as 
illustrated in Figure 4.17. For the FROG study a roll rate of 20°/sec is selected. Initially 
a lower roll rate of 5°/sec was used. However, the lower rate produced small values of 
both side force and yawing moment. A higher roll rate prevented machine resolution 
from becoming a factor in data reduction. The User Guide in Appendix A provides a 
more detailed description of the techniques used to gather the dynamic lateral-directional 
derivatives. 

The roll damping stability terms are presented in Table 4.10. Traditional 
design calculations from Roskam [Ref 12] were completed for comparison to CMARC 
and flight-test results. In addition, C-172 data is also presented. CMARC roll damping, 
Cip, is 50% higher than classical results, but right in line with parameter estimation. 
CMARC adverse yaw due to roll rate, Cnp, is approximately one fourth that obtained from 
either classical calculations or parameter estimation. Of note, CMARC does capture side 
force due to roll rate. This term is difficult to obtain from classical techniques and is 
usually assumed to be negligible. 


METHOD 

CONFIGURATION^ 

ROLL RATE TERMS 

cyp 

(per raaj 

C|p 

(per rao) 

^np 

(per rao) 

CMARC 

Panel Code 

Wing/Fuselage + HorzA/ert Stabs 

0.05 

-0.452 

-0.022 

Classical 

Design^ 

Wing/FuselageA/ert Tail 

0 

-0.300 

-0.072 

Parameter 

Estimation^ 

Flying Aircraft 

0 

-0.448 

-0.108 


C-172'' 

-0.037 

-0.47 

-0.03 


NOTES: 1) CGx=34.5% M.A.C. / 00^=8.6" from bottom of fuselage. 

2) Classical design calculations by Roskam's methods, after Ref. [12]. 

3) Parameter estimation from flight test data by Engdahl. 

4) 0-172 data from Ref. [12], 

Table 4.10 Comparison of FROG UAV Roll Rate Stability Derivatives. 
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The yaw rate-damping terms are presented in Table 4.11. Traditional 
design calculations from Roskam [Ref. 12] are presented for comparison to CM ARC and 
flight-test results. C-172 data is also listed. CMARC does an excellent job of predicting 
yaw damping. Yaw damping, Cnr, directly matches the result from parameter estimation. 
CMARC underestimates roll moment due to yaw rate, Cir, by 25 to 40 percent. And, side 
force due to yaw rate, Cvr, is overestimated by a factor of three compared to the 
parameter estimation model. No explanation is readily available to explain these errors. 


METHOD 

CONFIGURATION' 

YAW RATE TERMS 

Cvr 

(per rao) 

C,r 

(per rao) 

Cnr 

(per rao) 

CMARC 

Panel Code 

Wing/Fuselage + HorzA/ert Stabs 

0.337 

0.121 

-0.121 

Classical 

Design^ 

Wing/FuselageA/ert Tail 

0.140 

0.168 

-0.076 

Parameter 

Estimation^ 

Flying Aircraft 

0.110 

0.208 

-0.121 


C-172^ 

0.210 

0.096 

-0.099 


NOTES: 1) CGx=34.5% M.A.C. / CG 2 = 8 . 6 " from bottom of fuselage. 

2 ) Classical design calculations by Roskam's methods, after Ref. [ 12 ]. 

3) Parameter estimation from flight test data by Engdahl. 

4) C-172 data from Ref [12]. 

Table 4.11 Comparison of FROG UAV Yaw Rate Stability Derivatives. 

c. Lateral-Directional Control-Power Derivatives 


Aileron and rudder control-power derivatives are presented in Tables 4.12 
and 4.13. Classical design calculations are performed using equations 4.25 through 4.30 
from Roskam [Ref 15]: 





4.25 
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4.26 


C =C 

^'Sr '^rSr 


{Z^ cosa-l^ sin a) 


C, 


"A- 


^ (/^ COS a-Z^ sin or) 


4.27 


Cy^ « small 


4.28 


c„^=k*c,*q 


^Sa 


4.29 




& 


*g 

2;t 


^C' ^ 
y 


4.30 


where Cy 5 f, as, k and (C’isa/K) are empirical values from charts in Roskam [Ref. 15]. 
Additionally, flight-test parameter analysis and C-172 data are provided for comparison. 


METHOD 

CONFIGURATION^ 

AILERO 

CONTROL POWER 

Cysa 
(per rao) 

C|8a 

(perrao) 

^n6a 
(per rao) 

CMARC 

Panel Code 

Blended Wing-Fuselage/HorzA/ert Tails 

-0.021 

0.194 

-0.0121 

Classical 

Design^ 

Wing/FuselageA/ert Tail 

0 

0.213 

-0.0236 

Parameter 

Estimation^ 

Flying Aircraft 

0 

0.239 

-0.0261 


C-172^ 

0 

0.178 

-0.053 


NOTES: 1) CG„=34.5% M.A.C. / CG 2 = 8 . 6 ” from bottom of fuselage. 

2) Classical design calculations by Roskam's methods, after Ref. [12]. 

3) Parameter estimation from flight test data by Engdahl. 

4) C-172 data from Ref. [12]. 

Table 4.12 Comparison of FROG UAV Aileron Control-Power Derivatives. 
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In general, CMARC provides aileron control-power estimates in the 
correct direction and same order of magnitude as both the classical and parameter 
estimation techniques. Aileron roll control-power shows excellent correlation to both 
methods. However, CMARC under-estimates adverse yaw due to aileron deflection, 
Cn 5 a, by approximately 50%. Yawing moment due to roll rate, Cnp, is also 
underestimated by a similar margin. An investigation into the source of CMARC 
inaccuracies in the modeling of adverse yaw due to aileron deflection and roll rate is 
warranted. 

The rudder control-power derivatives obtained with CMARC show good 
correlation to classical design techniques. Side force due to rudder deflection, C\st, is 
15% greater than classical calculations but closely matches flight-test results. For both 
Cn6r and Ci6r, CMARC estimates fall between classical calculations and parameter 
estimation results. The stronger rudder control-power observed in the parameter 
estimation model is most likely due to the capturing of prop wash effects. 


METHOD 

CONFIGURATION^ 

RUDDER CONTROL POWER 

Cysr 
(per rao) 

C|5r 

(per rao) 

CnSr 
(per raa) 

CMARC 

Panel Code 

Blended Wing-Fuselage/HorzA/ert Tails 

0.0928 

0.0040 

-0.0453 

Classical 

Design^ 

Wing/FuselageA/ert Tail 

0.081 

0.0056 

-0.0341 

Parameter 

Estimation^ 

Flying Aircraft 

0.093 

0.0004 

-0.0785 


C-172^ 

0.187 

0.015 

-0.0657 


NOTES; 1) CGx=34.5% M.A.C. / 00^=8.6" from bottom of fuselage. 

2) Classical design calculations by Roskam's methods, after Ref. [12]. 

3) Parameter estimation from flight test data by Engdahl. 

4) C-172 data from Ref [12]. 

Table 4.13 Comparison of FROG UAV Rudder Control-Power Derivatives. 
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3. 


Summary' of CMARC Stability Derivative Analysis 


Table 4.14 lists the complete FROG aerodynamic model obtained from CMARC 
for the trim condition. In summary, CMARC produces reasonably accurate stability 
derivatives for an initial aerodynamic model. The greatest difficulty is encountered 
modeling side force due to sideslip and roll rate. The potential flow solution from 
CMARC fails to adequately capture the side force on the fuselage. Improved fidelity 
might be obtained by modeling fuselage flow separation with wake separation lines. 
CM.ARC also underestimates yawing moment due to roll rate and aileron deflection by 
large margins. 


LONGITUDINAL 

LATERAL-DIRECTIONAL 

Derivative 

Value 

Derivative 

Value 

Cl 

0.4295 

Cyp' 

-0.2493 

Cd 

0.065 

C|3 

-0.0630 

Cta 

4.845 

Cn3 

0.6300 

^Da 

0.2664 

^Yp 

0.0488 

^Ma 

-0.4126 

C,p 

-0.4514 

^Ladot 

1.420 

^np 

-0.0220 

^Madot 

-6.264 

Cvr 

0.3370 

Clg 

6.862 

C,r 

0.1210 

^Dq 

0 

Cnr 

-0.1210 

^Mq 

-11.78 

^Y5r 

0.0928 

Clsd 

0.4378 

Cl5r 

0.0040 

^D6e 

0.0092 


-0.0453 

^M6e 

-1.199 

^Yea 

-0.0206 


C^I5a 

0.1943 

^n6a 

-0.0121 


Table 4.14 Summary of CMARC Stability Derivatives for the NFS FROG UAV, 
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F. 


COMPARISON OF DYNAMIC AERODYNAMIC MODELS 


This section will discuss the dynamic response of the FROG UAV using the 
CMARC generated aerodynamic model. Specifically, modal frequency and damping are 
obtained through eigenvalue analysis. Elevator, aileron and rudder control response is 
found using a linearized state-equation model. Results are compared to the classical 
design and parameter estimation models. 

Dimensional stability derivatives are placed into linearized 4x4 models as 
outlined by Schmidt in Ref [14]. The short-period and long-period modes are obtained 
using the MATLAB eigenvalue decomposition routine. In a similar manner, the lateral- 
directional plant matrix yields the roll, spiral and Dutch-roll modes. Dynamic response 
due to step and doublet control inputs is obtained by assembling the complete linear 
system. This paper uses control deflection sign convention consistent with Figure 4.9 
from Ref [14]. The MATLAB "Isim" command produces the time based output of the 
linear system. MATLAB is used to plot the dynamic response. 


1. Longitudinal Dynamics 


a. Longitudinal Dynamic Modes 

The longitudinal response of an aircraft can be reduced to a series of four, 
first-order differential equations. They are typically written in matrix form based on the 
state variables uA^, a, q and 0. The linearized state-equations, in matrix form as 
developed by Schmidt in Ref [14], are listed in equations 4.31 to 4.37: 

{4=[.4]{x)+{sK 4.31 

{x]=[ulV a q ej 4.32 


[AAi.VM 
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4.34 




0 

-gcos ©0 




-gsin©o 

VM^ 


M, 

0 

0 

0 

1 

0 


4.35 


{B„]=[X^ or 


4.36 



V 

0 

0 

0 


0 0 0 
(F-Zj 0 0 
-M, 1 0 

0 0 1 


Inertial matrix 


4.37 


The non-dimensional stability derivatives generated by CMARC are 
converted to dimensional derivatives through the transformations listed in Table 4.15. 
Dimensional derivatives are more convenient because they will lead to time histories 
being expressed in seconds and frequency in rad/s. The MATLAB script "froguav.m" in 
Appendix D is used to automate these calculations. 

The linearized, longitudinal, state-equation is given by equation 4.31. 
Note that the plant and control matrices are obtained by pre-multiplying by [In] *• The 
eigenvalues of the plant matrix, [A], yield the longitudinal dynamic modes. The 
eigenvalues are obtained from the MATLAB script "froguav.m" listed in Appendix D. 
The results are presented in Table 4.16. 
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Term 

Description 

Units 

Xu 

mV dM 

s-1 

X. 

QS_ aco- 

da * 

ft-s~^ 

X . 

QS / c \ (9 Cd 

ft-s~^ 


m \2V/ (9(dc/2V) 

X 

Q S / c \ (9 Cd 

ft-s“^ 


m \2y) a(qc/2V) 

x. 

QS aco 

ft-s“^ 


m dS 

Zu 


s-i 

Za 

m aa 

ft-s“^ 

7 . 

_ Q S / c \ acL 

ft-s~^ 

"a 

m V2V/ a(ac/2V) 

7 

^ / c \ acL 

ft-s“ ^ 


m I 2 V/ a(qc/2V) 

Zi 

QS acL 

ft-s“^ 


m dS 

Mu 

Q S c dCm 
lyV aM 

ft-s“^ 

M 

QScaCm 

_ •> 

IVl^ 

ly da 

s 


QSC/ c \ aCm 


ly V2V/' a(dc/2V) 

s 

Mq 

QSC / C X aCn, 

c-1 

ly V2vJ a(qc/2V) 

s 

Mi 

QSc aCm 

Iv dS 

S 


Table 4.15 Dimensional Longitudinal Stability Derivatives from Ref. [14]. 
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METHOD 

configuration’ 

LONGITUDINAL DYNAMIC MODES 

Short- 

Period 

Phugoid 

COn 

Uao/s) 

C 

COn 

(raa/s) 

c 

CMARC 

Panel Code 

Wing/Fuselage/Horiz Tail 

5.15 

0.909 

0.358 

0.136 

Classical 

Design^ 

Wing/Horiz Tail - 58/5a=0.40 

5.90 

0.734 

0.407 

0.110 

Parameter 

Estimation^ 

Flying Aircraft 

4.67 

0.770 

0.397 

0.101 


C-172 

6.027 

0.685 

0.181 

0.116 


NOTES: 1) CGx=34.5% M.A.C. / CG2=8.6" from bottom of fuselage. 

2) Classical design calculations, after Ref. [12], 

3) Classical model from Papageorgio [Ref. 2] as modified through 
parameter estimation from flight test data by Engdahl. 

4) C-172 data from Ref. [12]. 

Table 4.16 Comparison of FROG UAV Dynamic Longitudinal Modes. 

Table 4.16 summarizes the longitudinal dynamic modes. The CMARC 
aerodynamic model provides a better estimate of short-period frequency than the classical 
design technique. The CMARC short-period frequency is 10% faster than that observed 
from parameter analysis and damping is 21% higher. The high damping in all three 
models will produce nearly deadbeat results. The CMARC phugoid, or long period, is 
10% slower than observed from flight-test with 40% more damping. Although, it must 
be pointed out that it is difficult to accurately capture the phugoid mode in flight-test. 
This is due to difficulty the external pilot has in accurately maintaining wings level trim 
during the extended period required for capturing phugoid data. Overall, the CMARC 
aerodynamic model produces satisfactory short-period and phugoid modal data for the 
development of closed-loop flight controls. Of note, CMARC provides a more accurate 
prediction of the important short-period natural frequency. The larger error in phugoid 
frequency is less important to the design of a closed loop controller. 
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b. Longitudinal Dynamic Response to Control Input 

Following dynamic mode analysis, the response to an elevator step-input 
and a doublet-input is modeled using MATLAB. Dynamic response for the CMARC, 
classical and parameter analysis models are overlaid for comparison. The linear system 
is set up using the MATLAB script "froguav.m” in Appendix D. The "Isim" command 
outputs the time-based response from a time-based control input vector. A -2° (TEU) 
elevator step-input is selected to keep the response in the linear region. A 5° elevator 
doublet with a 1.2 second period is used to excite the short-period mode. The selected 
doublet period closely matches the short-period mode. Dynamic response is displayed in 
Figures 4.18 and 4.19. 

The FROG UAV response to a -2° (TEU) step elevator input is displayed 
in Figiire 4.18. All three aerodynamic models are presented for comparison. The model 
adjusted by parameter analysis shows a much larger response than either the CMARC or 
classical design models. The CMARC model produces a final pitch angle change after 4 
seconds of 35° versus 45° for the flight-tested model. On closer review, it is noted the 
CmSe =-1.199 from CMARC is considerably less than the CmSe = -1.621 from flight-test. 
It should be noted that the parameter estimation model compensates for real world factors 
including; air load distortions, sensor measurement errors and prop wash. In other words, 
the CMARC model assumes uniform displacement and perfect sensors while the flight- 
test model is empirically fit to a measured response. 

The response to a 5° elevator doublet is displayed in Figure 4.19. All 
three models show a similar frequency and a high degree of damping. As expected, the 
parameter estimation model produces a larger response. The magnitude of the CMARC 
response is approximately two thirds of that observed from the flight-test model. Again, 
this is due to the much larger value of Cm 5 e obtained from empirically fitting the 
parameter estimation model to observed aircraft response. Modeling prop wash with 
CMARC should produce improved results. 

In summary, the CMARC aerodynamic model showed satisfactory 
longitudinal control response in frequency and damping. It is recommended that CmSe be 
adjusted to increase elevator response. 


65 


4 


LONGITUDINAL RESPONSE TO A 2 DEGREE ELEVATOR STEP INPUT 
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Figure 4.18 FROG UAV Dynamic Response to a - 2 ° (TEU) Elevator Step Input 
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LONGITUDINAL RESPONSE TO A 5 DEGREE ELEVATOR DOUBLET 






Figure 4.19 FROG UAV Dynamic Response to a 5® Elevator Doublet 
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2 . 


Lateral-Directional Dynamics 


a. Lateral-Directional Dynamic Modes 

The lateral-directional response of an aircraft can also be reduced to a 
series of four first-order differential equations. They are typically written in matrix form 
based on the state-variables P, p, (j), and r. The linearized state-equations in matrix form 
as developed by Schmidt in Ref [14] are listed in equations 4.38 to 4.44: 




{A=\p p <!’ T 
M=[/.]-'kl 

{5}=[/,]-'k} 



0 


K 

1 


geos 00 (K + v) 
0 4 

0 0 

0 




0 

0 



V 

0 

0 

0 


0 0 0 
(f-zJ 0 0 

-M, 1 0 

0 0 1 


Inertial matrix 


4.38 

4.39 

4.40 

4.41 

4.42 

4.43 


4.44 
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As with the longitudinal axis, the stability derivatives generated by 
CMARC are dimensionless. However, it is more convenient to use dimensional 
derivatives. They will lead to time histories being expressed in seconds and frequency in 
rad/s. The lateral-directional stability derivatives are turned into dimensional derivatives 
through the transformations listed in Table 4.17 from Ref [14]. The MATLAB script 
"froguav.m" in Appendix D is used to automate these calculations. 



Table 4.17 Dimeosional Lateral-Directional Stability Derivatives from Ref. [14] 
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The linearized, lateral-directional state-equation is given by equation 4.38. 
Note that the plant and control matrices are obtained by pre-multiplying by [In] '- The 
dimensional derivatives from Table 4.16 populate the linearized 4x4 plant matrix and 
control matrices. The plant matrix, [A], is used to obtain the lateral-directional dynamic 
modes. The eigenvalues are obtained from MATLAB using the "froguav.m" script in 
Appendix D. The results are presented in Table 4.18. 


METHOD 

CONFIGURATION^ 

LATERAL-DIRECTIONAL DYNAMIC MODES 

Dutc 

h-Roll 

Roll 

Spiral 

COn 

(rao/s) 

c 

COn 

(rao/s) 

(On 

(raa/sj 

CMARC 

Panel Code 

Wing/Fuselage Group 
plus HorizA/ert Stab Group 

2.54 

0.13 

-3.80 

0.000 

Classical 

Design^ 

Wing/Fuselage/Horiz Tail 

2.48 

0.09 

-2.83 

0.065 

Parameter 

Estimation^ 

Flying Aircraft 

4.27 

0.14 

-3.97 

0.090 


C-172^ 

3.38 

0.20 

12.43 

0.01 


NOTES: 1) CGx=34.5% M.A.C. / CGz=8.6" from bottom of fuselage. 

2) Classical design calculations after Ref. [12], 

3) Classical model from Papageorgio [Ref. 2] as modified through 
parameter estimation from flight test data by Engdahl. 

4) C172 data from Ref. [12], 

Table 4.18 Comparison of FROG UAV Lateral-Directional Dynamic Modes. 

The lateral-directional dynamic modes are summarized in Table 4.18. The 
CMARC aerodynamic model provides significantly better estimates of frequency and 
damping than the classical design technique for both the Dutch-roll and roll modes. The 
CMARC Dutch-roll natural frequency, although better than the classical response, is 40% 
lower than predicted by parameter estimation. Adjustments should be made to the 
CMARC aerodynamic model to improve Dutch-roll response. 
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The primary contributor to the Dutch-roll frequency is weathercock 
stability, Cnp. A review of Table 4.9 shows that Cnp obtained from CMARC is 
considerably less than the value obtained through parameter estimation. One potential 
source of error is the lack of modeling the increased dynamic pressure due to prop wash. 
The FROG model should be re-worked to include a prop disk to investigate the ability of 
CMARC to capture prop wash effects. 

The CMARC model predicts a neutral spiral mode. This is most likely 
due to a combination of equal and opposite Cir and Cnr ratios. When the Cir/Cnr ratio is 
changed to -2 (approximately that of the analytical and parameter estimation models), the 
spiral mode comes out close to the analytical model. The weak roll due to yaw-rate, Cir, 
in the CMARC solution seems to be the main source of error in the spiral mode. The 
value will need to be modified to produce an acceptable spiral response. 

b. Lateral-Directional Dynamic Response to Control Input 

Following dynamic-mode analysis, the response to aileron and rudder 
inputs is modeled using MATLAB. Dynamic response for the CMARC, classical and 
parameter-analysis models are overlaid for comparison. The linear system is set up using 
the MATLAB script "froguav.m" in Appendix D. The "Isim" command outputs the time- 
based response from a time-based input control vector. 

The lateral response to a +2° (right wing down) step aileron deflection is 
shown in Figure 4.20. CMARC steady-state roll-rate is less than that predicted by the 
other two models. This is expected considering the CMARC model has the lowest ratio 
of aileron control-power to roll-damping. In addition, the side-force terms, Cyp and Cysa, 
modeled by CMARC both fight roll-rate. They are not included in the other models. 
Perhaps they should be assumed small and set to zero. CMARC does a significantly 
better job of modeling the Dutch-roll response. The roll-rate and sideslip traces from the 
classical design model clearly show excess excitation of the Dutch-roll mode. The 
CMARC model shows a better correlation to the parameter estimation model for Dutch- 
roll amplitude and damping. As predicted by the eigenvalues. Dutch-roll frequency is 
slower than the parameter estimation model. 

The lateral-directional response to a 5° aileron doublet is shown in Figure 
4.21. The 1.5 second period clearly excites the Dutch-roll mode in all three models. As 
with the step-input response, the CMARC model provides a significantly better match to 
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Dutch-roll excitation than the classical design model. Amplitude and damping show a 
close correlation to the parameter estimation model. A lower Dutch-roll frequency is 
evident in the CMARC sideslip trace. 

The lateral-directional response to a +2° (nose left) rudder deflection is 
shown in Figure 4.22. The CMARC model demonstrates a close match to the classical 
design model for Dutch-roll frequency and damping. CMARC shows higher sideslip 
excitation than the classical design model. A slower Dutch-roll frequency is evident in 
both the CMARC and classical design models. 

The lateral-directional response to a 5° rudder doublet deflection is shown 
in Figure 4.23. Again, the 1.5 second doublet excites the Dutch-roll mode in all three 
models. CMARC provides a similar response to the classical design model. As 
expected, Dutch-roll natural frequency is about 40% slower than the parameter estimation 
model with similar damping. 

In summary, the CMARC lateral-directional model provides FROG 
dynamics that are similar to the classical design calculations. However, current CMARC 
lateral-directional dynamics are not adequate for closed-loop controller or autopilot 
design. Minor adjustments will be required to better match observed flight 
characteristics. The lateral-directional model should be modified to include a propeller 
disk. Higher dynamic pressure over the tail surfaces should provide a stronger Cnp 
derivative and result in faster Dutch-roll frequency. In addition, the FROG model should 
be adjusted to provide a higher aileron control-power to roll-damping ratio. This 
modification will improve roll-rate response. 
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Figure 4.20 FROG UAV Dynamic Response to a +2° (RWD) Aileron Step Input 
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LATERAL-DIRECTIONAL RESPONSE TO A 5 DEGREE AILERON DOUBLET 






Figure 4.21 FROG UAV Dynamic Response to a 5° Aileron Doublet. 


74 





































4 


LATERAL-DIRECTIONAL RESPONSE TO A 2 DEGREE RUDDER STEP INPUT 
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Figure 4.22 FROG UAV Dynamic Response to a +2° (TEL) Rudder Step Input 


75 


































lateral-directional response to a 5 DEGREE RUDDER DOUBLET 




-!- T -----,- 1 - 




- CMARC 


- 



Classical 

- 

1 

_I_ \ _1_ 

Parameter Est 



012345678 





Figure 4.23 FROG UAV Dynamic Response to a 5° Rudder Doublet 
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V. CONCLUSIONS AND RECOMMENDATIONS 


CMARC is a DOS personal-computer hosted panel-code adopted from the NASA 
Ames PMARC code. AeroLogic, Inc., created CMARC by converting PMARC 
FORTRAN 77 source code into the C language. Significant memory management and 
command-line enhancements were also added. CMARC solves for inviscid, 
incompressible flow over complex three-dimensional bodies. Emphasis in this study is 
placed on expanding the use of CMARC to develop a full aerodynamic model of the 
Naval Postgraduate School FROG UAV. The CMARC stability derivatives are 
compared to models derived from classical design calculations and parameter estimation. 
In addition, pitot-static and angle-of-attack sensor position corrections are obtained 
through CMARC analysis. 

The LOFTSMAN and POSTMARC portions of the Personal Simulation Works 
software suite are used exclusively for the pre-process modeling and post-process 
visualization of CMARC files. The LOFTSMAN capability to automatically format and 
generate CMARC input patches is an enhancing characteristic. Functionality should be 
added to allow the modeling of wing tip ribs that are not parallel to the aircraft centerline. 

POSTMARC is an excellent tool for visualizing CMARC output files. The 
capability to create streamlines and perform boundary-layer calculations external to 
CMARC is extremely useful. However, much time could be saved if POSTMARC 
maintained previous settings and selections following translations, rotations and re¬ 
scaling. Additionally, a capability to overlay multiple data types is desired. 

CMARC off-body flow-field analysis is useful for both static source and angle-of- 
attack sensor position corrections. In-flight measurements may be corrected using look¬ 
up tables or through curve-fits of CMARC-derived data. Flight testing is recommended 
for validation of sensor corrections obtained from the CMARC off-body analysis. 

For the static longitudinal analysis, CMARC produces accurate values for 
and C^a ^ slightly high value of Clo- Elevator control-power (Cmsc) from CMARC is 
considerably weaker than the value obtained from parameter estimation. CMARC- 
derived pitch-damping is stronger than the pitch-damping from both classical design 
calculations and parameter estimation. 

Longitudinal dynamic-mode analysis shows an acceptable match for short-period 
and long-period frequency and damping. As expected, a lower dynamic response to 
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elevator control inputs is observed when compared to the parameter estimation model. 
However, in all cases, the CMARC aerodynamic model demonstrates better dynamic 
response than the classical design model. A propeller disk should be added to the 
CMARC model in an attempt to capture prop-wash effects on the horizontal stabilizer 
and elevator. Increased elevator control-power would provide a better match with 
observed flight characteristics. 

In general, CMARC provides satisfactory lateral-directional derivatives. Roll¬ 
damping closely matches the value obtained through parameter estimation. The greatest 
difficulty encountered is in modeling side-force due to both sideslip and roll-rate. The 
potential flow solution from CMARC fails to adequately capture the side force on a 
slender-body fuselage. Additional wakes, placed along separation lines, may improve 
fuselage side-force prediction. 

Dynamic response to rudder control input shows close correlation in amplitude 
and damping to the parameter estimation model. However, both the CMARC and 
classical models of the FROG UAV produced a slower Dutch roll frequency. Steady 
state roll-rate obtained from the CMARC model is somewhat slower than either the 
classical design or parameter estimation models. 

Overall, the CMARC panel code is found to be suitable for aerodynamic 
modeling of the Naval Postgraduate School FROG UAV. CMARC-derived stability 
derivatives are sufficiently accurate for incorporation into an initial aerodynamic model. 
The CMARC aerodynamic model demonstrated better longitudinal dynamic response 
than the classical design model. Lateral-directional response is similar to that obtained 
from classical design techniques. Adjustment through comparison with flight-test data is 
still required to optimize the CMARC model. Future studies should concentrate on 
improving CMARC modeling of fuselage side force through the addition of separation 
wake lines. Additionally, the propeller disk should be modeled in an attempt to capture 
the effects of increased dynamic pressure over the horizontal and vertical tail surfaces. 
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A. 


INTRODUCTION 


1. CMARC Panel Code and Limitations 

CMARC is a PC hosted panel code. The code is useful for gathering aircraft 
stability data subject to the limitation of inviscid, potential flow solutions. Panel codes 
assume fully attached flow. Therefore, care must be taken to gather stability data in the 
linear region of the lift curve slope where separation effects are minimized. The results 
obtained from CMARC should always be spot checked with classical design calculations. 

CMARC and PMARC have been shown to produce functionally equivalent 
results. The guide describes the use of DOS batch files for automating CMARC data 
runs on a PC. PMARC, hosted on an SGI machine, will require the use of UNIX batch 
conunands. 

2. Overview of Obtaining Stability Derivatives with CMARC 

An overview of the methods to obtain stability derivatives follows. The detailed 
description and data reduction methods can be found in the appropriate section of the user 
guide. 

For the static stability derivatives, the CMARC model is run at two different 
angles-of-attack and one sideslip angle. The static solutions are obtained without control 
surface deflections. The slopes of the force and moment coefficients are then taken to 
produce the C^j and Cn,„ longitudinal derivatives and the Cyp, C|p and Cnp lateral- 
directional derivatives. 

For the control power derivatives, the model is run at the trim condition with 
successive elevator, aileron and rudder control deflections. The difference between the 
results with and without control surface deflections yield the Clscj C^se and Cps^ 
longitudinal and the Cysr, C, 5 r, C^j and C„ 5 a lateral-directional control power 
derivatives. 

Development of the damping derivatives is not as straight forward. For the 
dynamic derivatives, motion is enabled around the center of gravity. With motion turned 
on, it is extremely important to ensure that the global origin and C.G. are co-located. If 
they are not, a rotation around a specific axis will also create a translation at the C.G. 
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For the longitudinal damping derivatives, the model is run with oscillating vertical 
plunging motion to obtain the Cl and a-dot terms. The lift and pitching moment 
coefficients are broken into real (in phase with AOA) and imaginary (out of phase with 
AOA) components. The imaginary components are due to a-dot effects. Next, the model 
is run with oscillating pitch motion to obtain the combined a-dot and pitch rate terms. 
Subtracting the a-dot influence obtained fi-om the plunging motion isolates the pitch rate 
damping term from the pitch motion. 

For the lateral-directional analysis, the P-dot terms are generally negligible. This 
allows the model to be run with just oscillating roll and yaw motion. As with the 
longitudinal test case, the imaginary or out of phase component yields the combined p- 
dot and rate damping data. With the p-dot terms assumed negligible, the oscillating 
motion yields the Cy, C„ and C„ roll and yaw rate terms directly. 

3. Coefflcient Normalization and Stability Axis System 

CMARC contains built-in functionality to integrate forces and moments over the 
surface of a body. Forces and moments are automatically normalized into non- 
dimensional coefficients based on the mean aerodynamic chord, reference wing area, 
semi-span and center of gravity location in the CMARC BINP9 input line. CMARC 
outputs coefficients in both the wind and body axes. Of note, CMARC uses the semi¬ 
span to normalize rolling and yawing moment coefficients. Most texts on stability and 
control, including Roskam's "Aircraft Flight Dynamics" and Etkin's "Dynamics of 
Flight," normalize rolling and yawing moments by span. This user guide will normalize 
roll and yaw moments by span. Table A.l summarizes the factors for normalizing 
moments and angular rates. All rolling and yawing moment coefficients presented in this 
guide have been normalized with span by dividing the CMARC output by a factor of two. 
Table A.l also indicates the characteristic time, t*, employed for angle rate data 
reduction. 

In addition to differences in normalizing moments, CMARC uses the typical CFD 
axes system shown in Figure A.l. This user guide will perform all calculations in the 
stability axes system as illustrated in Figure A.l. The sign of CMARC roll and yaw 
moments need to be reversed. The direction for positive control deflections is also shown 
in Figure A.l. All control surfaces are patched with positive deflections using the 
convention in Figure A.l. 
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NORMALIZING 


CHARACTERISTIC 

MOMENTS 

PARAMETER’ 

RATES 

TIME 

II 

^1 

b 

II 

r= ¥ 

M = C.qSc 

c 

^ rc 

r = - 

2u^ 


N = C^qSb 

b 

II 

¥ 


Note: 1) CMARC normalizes roll and yaw coefficients with b/2. 


Table A.l Normalized Moment and Rate Equations. 

B. MODEL VALIDATION 

1. Initial Solution at the Trim Condition 

Once a CMARC model is built up, a reference solution is obtained at the desired 
trim condition. Initial hand calculations or flight test data can be used to select the 
estimated cruise angle-of-attack. Use POSTMARC to visualize the reference solution. 
Check to make sure the pressure distribution is consistent with expected results. 
Stagnation zones (Cp=l) should be evident on the leading and trailing edges of the flight 
surfaces and aerodynamic bodies. Lower pressure should be observed on the top of 
lifting surfaces and tip loss should be observed near the wing tips. Wakes should trail 
from the expected surfaces. 

If abnormal pressure distribution is observed, check for inverted panel orientation. 
This is accomplished by changing the color for positive and negative panel orientations 
using the "View » Color" pull-down menu. The default color is white for both 
orientations. Change the inside to a contrasting color such as red. All panels should face 
in the positive out direction. If a patch has panels oriented inward, change the orientation 
of the patch with the "IREV" term in the initial patch definition. 

Check that the model is symmetric around the lateral and directional axes with 
sideslip set to zero. Residual rolling or yawing moment without a control deflection is an 
indication of an asymmetric model. 
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CMARC Axis System 



Stability Axis System 


Figure A.l CMARC CFD Axis System Compared To Stability Axis System. From 
References (3) and [14] Respectively. 


86 













Pitching moment should be checked for a reasonable value. Don't forget to verify 
that the origin of the global coordinate system is set to the center-of-gravity (C.G.). It is 
important to remember that the damping derivatives will be obtained by inducing motion 
around the center-of-gravity. If the C.G. and origin are not collocated, angular motion 
will generate undesired translation at the C.G. 

Next, the lift coefficient from the initial solution is checked to make sure it is in 
the "ballpark." If a major discrepancy is noted, the input file will have to be checked 
closely for error. Be sure to check the normalizing factors such as wing area, reference 
chord and semi-span. If every thing looks good, try increasing or decreasing the angle- 
of-attack to see if a slight adjustment will bring the hft coefficient in line with the hand 
calculation. Once the model exhibits symmetry with a reasonable lift distribution, you're 
ready to start gathering data. 
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c. 


LONGITUDINAL STABILITY DERIVATIVES 


This section will describe CMARC methods for the development of longitudinal 
stability derivatives to include the static, rate damping and control power derivatives. As 
stated earlier, the results obtained from CMARC should be checked against classical 
design calculations. The potential flow analysis performed by CMARC does not provide 
accurate viscous drag values. Total drag can be estimated from the flight test power-off 
glide ratio or cruise thrust required. If this data is unavailable, then empirical data from 
publications such as Hoemer's "Fluid Dynamic Drag" can be used to estimate total drag. 

For the longitudinal analysis, only half the model is analyzed. This cuts the 
model size in half, resulting in much quicker solutions. The symmetric calculation mode 
is selected by setting both RSYM=0.0 and IPATSYM=0 in the CMARC input file. 
Remove the vertical tail patch if it interferes with wing or fuselage wakes. Figure A.2 
shows the FROG model configuration used to find the longitudinal derivatives. A rigid 
wake is selected that is continuous from wingtip to wingtip. As will be seen in a later 
section, the vertical tail is added for the lateral-directional solutions. This requires the use 
of a modified wing wake to avoid the vertical tail. 




Blended Wing and Fuselage with Horizontal Tail 


Figure A.2 FROG UAV Model for Obtaining Longitudinal Stability Derivatives. 
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1. 


Static Longitudinal Stability Derivatives 


a. Static Longitudinal Stability Equations 


Three basic longitudinal stability derivatives can be measured with just 
two runs of the CMARC model. The model is first analyzed at an angle-of-attack 
corresponding to the estimated trim condition. For the FROG case, a=0° is selected for 
the first run based on hand calculations. A second CMARC run is conducted with angle- 
of-attack incremented by one or two degrees. C^, and Cn, are then extracted manually 
from the data files. The slope of Cl, C^ and C^, versus angle-of-attack provides the Clo, 
Cdo and Cn,„ longitudinal derivatives. Additional solutions should be obtained to check 
for lift slope linearity. In addition, is calculated fi-om the lift curve slope and trim lift 
coefficient. Equations A.l through A.4 are used for these calculations. 


Cl — 


-Cij J 180 


= 


C = 

'^m.r 


(«2 - 

■«l) 

n 

{Cd2 

-^Dl) 

,180 

(«2 

-«i) 

71 


-c„,J 

1,180 

(«2 

-«l) 

71 


per radian 


A.1 


A.2 


A.3 


a: 


trim 


= + 






*- degrees 

n 


A.4 


b. Sample Static Longitudinal Stability Data Reduction 

Table A.2 presents CMARC solutions for the FROG UAV at two angles- 
of-attack. Sample calculations for obtaining the static derivatives are demonstrated 
below. The calculations are easily implemented in a spreadsheet or with a MATLAB 
script. 
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RUN# 

AOA 

CL 

CD 

CY 

C m 

C n 

C 1 

1 

0° 

0.4300 

0.0277 

0.0000 

0.0444 

0.0000 

0.0000 

2 

2° 

0.5991 

0.0370 

0.0000 

0.0300 

0.0000 

0.0000 


Table A.2 FROG UAV Static Longitudinal Stability Data. 


Sample Calculations: 


C, 


_ (0.5991-0.4300)^180 




( 2 - 0 ) 


7t 


= 4.8444 per radian 


^ (0.0370-0.0277)^180 

Cn - 7 - X --= 0.2664 per radian 

(2-0) 7t ^ 


^ (0.0300-0.0444)^180 

Cm = - 7 -^- - -= -0.4125 per radian 

( 2 - 0 ) t: ^ 




trim 


^ ^ (0.4295-0.4300) . ^ _o.0059» » 0» 

4.8444 


7t 


A.5 


A.6 


A.7 


A.8 


c. Sample Total Drag Coefficient Calculation 

The final static longitudinal parameter required is total aircraft drag. Cp 
plays an important role in long period aircraft dynamics. Unfortunately, potential flow 
panel codes such as CMARC do not provide accurate total drag estimates. They can 
provide good induced drag predictions. In addition, if equipped with a boimdary layer 
code like that contained in CMARC, they can provide integrated skin fiiction results. 
However, a large total drag contribution in the form of separation drag is not accounted 
for. Total drag estimates are made below using the two simple techniques shown in 
Equations A.9 to A. 12. The first method is based on the flight-tested glide ratio. The 
second is based on cruise power required and estimated prop efficiency. Note that the 
selected prop efficiency is relatively low due to the FROG UAV small propeller diameter. 
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high RPM and pusher configuration. The two methods provide drag predictions within 
10% of each other. The results are averaged to Cd=0.065. 


Method 1: Lift-to-Drag Ratio (L/D=7 fi'om flight test) 


% = 7 = i = = lbs 

/D 111 


qS 


9.67 lbs 


0.5*0.002327 lb ■ ! ft* * ffls^*\1.51 ff 


= 0.0611 


A.9 


A.IO 


Method 2: Cruise Power Setting (HP=5, T|p=0.35) 

^ _ D _Tj, _ HP^*550 ftlbs/s/HP*qp/V 
qS~ qS~ qS 


A.11 




5 HP*550 ft-lbs/si HP* 035/SS ftls 
0.5*0.002327 lb ■ s^! ft* * ft^ls^*\1.51 ft^ 


0.069 


A. 12 


2. Longitudinal Control Power Stability Derivatives 


CL Longitudinal Control Power Equations 


The elevator control power derivatives are obtained by substituting a 0° 
defiection horizontal tail patch for one with positive elevator deflection. Only one run is 
required. For the FROG UAV study, +5° (TED) deflection is used. The difference 
between the trim condition and the deflected value is divided by the elevator deflection as 
shoAvn below. Note that Cds^ fi'om CMARC only includes induced drag due to elevator 
deflection: 


180 


^ _se 


n 


per rad 


= 


^DSe2 ^^Sei 180 


Sei -Sei 


n 


per rad 


A.13 


A.14 
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A.16 


5e° 


trim 


^^@trim{Se=(!i)) ^ 180 
n 


mse 


deg 


b. Sample Longitudinal Control Power Data Reduction 


Table A.3 presents CMARC solutions for two elevator deflections. 
Sample calculations for obtaining elevator control power are demonstrated below. The 
calculations are easily implemented in a spreadsheet or with a MATLAB script. 


C 


(0.4641-0.4259)^ 180 


j =- 7 -7-- ■-= 0.4378 per radian 

^5e (5-0) Tt ^ 


^ (0.0178-0.0170)^180 ^ 

C[)^ =- --r- -* -= 0.0092 per radian 

^ (5 “ 0 ) 71 


A.17 


A.18 


C 




_ (-0.0942-0.0104)^ 180 


(5-0) 


71 


= -1.199 per radian 


A. 19 




tnm 


(0.0104) ^ 180 
-1.199 71 


= -0.50° 


A.20 


RUN# 

5, 

CL 

CD 

CY 

C m 

C n 

C 1 

1 

0° 

0.4259 

0.0170 

0.0000 

0.0104 

0.0000 

0.0001 

2 

5° 

0.4641 

0.0178 

0.0000 

-0.0942 

0.0000 

0.0001 


Table A.3 FROG UAV Elevator Control Power Data. 

























3. Longitudinal a-dot Damping Derivatives 


a. Longitudinal a-dot Derivative Methods and Equations 

Two motions selected are to develop the longitudinal rate damping 
derivatives. A sinusoidal plimging motion in the z-axis is used to extract the a-dot terms. 
Then an oscillatory pitching motion is used to obtain the combined a-dot and pitch rate 
terms. The a-dot are then subtracted to yield the pitch rate damping. All motion for 
FROG data gathering is conducted at a frequency of In rad/s, which equates to a reduced 
frequency of k=0.0595 for this configuration and trim airspeed. 

The sinusoidal plunging motion illustrated in Figure A. 3 is used to isolate 
the a-dot derivatives. Z-axis plunging motion is controlled with the CMARC BINP8B 
input file line. A frequency of 27i rad/s and an amplitude of C/2=10 inches are selected 
for the FROG study. Greater amplitude was initially selected, but the induced angle of 
attack caused the rigid wake to impact the horizontal tail. An example BINP8B input line 
is shown below; 

&BINP8B DXMAX=0.0, DYMAX=0.0, DZMAX=10.0 

WTX=0.0, WTY=0.0, WTZ=6.283, &END 

The number of time steps and time step interval is chosen to create a nice 
sinusoidal output through at least two cycles of motion. In the FROG study, a plunging 
motion frequency of (o=2n rad/s or one cycle/sec is utilized. Fifty time steps are chosen 
with an interval of 0.05 seconds, which creates 2.5 cycles of plunging motion. After the 
solution is obtained, the "total coefficient" data is extracted for plotting. The data can be 
picked out manually, or a data retrieval program can be created for the task. For the 
FROG study, the data is extracted manually and pasted into an "Excel" spreadsheet for 
plotting. In addition to the CMARC output, plunging motion phase angle and induced 
angle-of-attack are calculated using Equations A.21 and A.22. Induced angle-of-attack is 
used to find the phase angle of the response with respect to angle-of-attack. Table A.4 
shows representative FROG data for 20 time steps. 
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^plunge timestep*^^ * ^plunge * ^ degrees 


^induced COS(^ plunge ) 


^plunge 




^plunge ^ degrees 


\\Tiere: 


<t>piunge - plunging motion phase angle 

u^induced • induced angle-of-attack from the plunging motion 

^timestep * time Step 

dt - time step interval 
©plunge - plunging frequency 
Apiunge - plunging amplitude 
Uq - reference free stream velocity 


A.21 


A.22 


Step 

< t > ( deg ) 

^Induced 

CL 

CD 

CY 

C m 

C n 

c l 

0 

0 

- 3.41 

0.0001 

0.0327 

0.0000 

0.0173 

0.0000 

0.0000 

1 

18 

- 3.24 

0.1318 

0.0300 

0.0000 

0.0744 

0.0000 

0.0000 

2 

36 

- 2.76 

0.1928 

0.0299 

0.0000 

0.0485 

0.0000 

0.0000 

3 

54 

- 2.00 

0.2631 

0.0316 

0.0000 

0.0294 

0.0000 

0.0000 

4 

72 

- 1.05 

0.3451 

0.0318 

0.0000 

0.0218 

0.0000 

0.0000 

5 

90 

0.00 

0.4336 

0.0297 

0.0000 

0.0169 

0.0000 

0.0000 

6 

108 

1.05 

0.5216 

0.0252 

0.0000 

0.0126 

0.0000 

0.0000 

7 

126 

2.00 

0.6009 

0.0189 

0.0000 

0.0087 

0.0000 

0.0000 

8 

144 

2.76 

0.6636 

0.0123 

0.0000 

0.0067 

0.0000 

0.0000 

9 

162 

3.24 

0.7029 

0.0069 

0.0000 

0.0089 

0.0000 

0.0000 

10 

180 

3.41 

0.7160 

0.0042 

0.0000 

0.0139 

0.0000 

0.0000 

11 

198 

3.24 

0.7009 

0.0048 

0.0000 

0.0232 

0.0000 

0.0000 

12 

216 

2.76 

0.6594 

0.0085 

0.0000 

0.0344 

0.0000 

0.0000 

13 

234 

2.00 

0.5952 

0.0142 

0.0000 

0.0474 

0.0000 

0.0000 

14 

252 

1.05 

0.5141 

0.0204 

0.0000 

0.0596 

0.0000 

0.0000 

15 

270 

0.00 

0.4254 

0.0255 

0.0000 

0.0704 

0.0000 

0.0000 

16 

288 

- 1.05 

0.3370 

0.0286 

0.0000 

0.0762 

0.0000 

0.0000 

17 

306 

- 2.00 

0.2582 

0.0297 

0.0000 

0.0776 

0.0000 

0.0000 

18 

324 

- 2.76 

0.1964 

0.0292 

0.0000 

0.0734 

0.0000 

0.0000 

19 

342 

- 3.24 

0.1561 

0.0283 

0.0000 

0.0713 

0.0000 

0.0000 

20 

360 

- 3.41 

0.1398 

0.0277 

0.0000 

0.0727 

0.0000 

0.0000 


Table A.4 Sample Plunging Motion Data for 20 Time Steps 
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Figure A.3 Vertical Plunging Motion for Obtaining the a-dot Derivatives. 

The a-dot stability derivatives are extracted using the methods outlined in 


Etkin [Ref. 13]. First, the Cl and Cn, responses are plotted as a function of plunging 
phase angle. Induced angle-of-attack is also plotted on the right hand vertical axis. 
Figure A.4 is a representative plot of CMARC data for the FROG UAV study. The 
response to plunging motion can be broken up into real and imaginary components. 
Figure A.5 is a graphical representation of the real and imaginary components. The out- 
of-phase (imaginary) portion of Cl or C^ is the a-dot contribution. It is normalized by 
dividing by the amplitude of and the reduced frequency. The phase angle is 

measured between the lift or pitching moment response and the induced angle-of-attack. 
Equations A.23 through A.27 are used to solve for CLadot 
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NPS FROG UAV 

Vertical Plunging Motion: coplunge=2pi rad/s, A=cbar/2 



Figure A.4 Representative Vertical Plunging Motion Data for the FROG UAV. 


Im 



Figure A.5 Complex Response to Vertical Plunging Motion. 
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Angle-of-Attack (degrees) 



















































The phase angle between the Cl (or C^) response and the induced angle- 


of-attack is measured graphically or through numerical curve fit techniques. If the 
graphical method is used, the graph scales may be narrowed aroimd the area of interest to 
facilitate the phase angle measurement. For instance. Figure A.4 shows that lift 
coefficient and angle-of attack are nearly in-phase. When the x-axis scale is zoomed into 
the peak of the curves, one can see that lift coefficient leads angle-of-attack by one 
degree. Alternatively, the phase angle may be measured fi'om where the parameter 
crosses the trim value. 



A.23 


C 



A.24 



o 


A.25 



o 


A.26 


A.27 


Where: (j) - Phase angle between Cl (or CJ response and 


^induced - Induced angle-of-attack from the plunging motion 

Clo - Amplitude of lift coefficient response 

Cmo - Amplitude of pitching moment coefficient response 

ttg - Amplitude of the induced angle-of-attack 

k - Reduced frequency 

Uq - Reference free stream velocity 


b. Sample Longitudinal a-dot Data Reduction 
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Sample data reduction is presented below in Equations A.28 through A.32. 
The phase angle between the response and the induced angle-of-attack was measured 
graphically by zooming the graph axes to expand the area of interest. 


q/ 20 m. 

k = r / 5 = 0.0595 

Uq 1056 in Is 


C, 


= Ci sm{<p)= -•^^^^^ *sin(l°)= 0.005028 


C„ 


= sin(<;5) = (0.0776 - 0.0087)* sin(- 140°)= -0.0222 


A lOinll _Ar>cn« 1 A 10 

cCq— ^plunge ~ ■ / 2;t — 0.0595 rad —3.41 

Un \056inls n 


C 




A^Lo\ 0.005028 


d a^k 0.0595*0.0595 


420 


C 


1 - 0-0222 


aok 0.0595*00595 


= -6.264 


A.28 

A.29 

A.30 

A.31 

A.32 
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4. 


LoDgitudinal Pitch Rate Dampiag Derivatives 


a. Longitudinal Pitch Rate Damping Methods and Equations 

A sinusoidal plunging motion in the z-axis was used to extract the a-dot 
terms. Now, an oscillatory pitching motion is used to obtain the combined a-dot and 
pitch rate terms. The a-dot influence is then subtracted to yield pitch rate damping. All 
motion for FROG data gathering is conducted at a frequency of In rad/s, which equates 
to a reduced frequency of k=0.0595 for this configuration and trim airspeed. 

The sinusoidal pitch rate motion illustrated in Figure A.6 is used to isolate 
the combined a-dot and pitch rate influence. Oscillating pitch motion is controlled with 
the CMARC BINP8A input file line. A frequency of 2tz rad/s and an amplitude of ±1.5° 
are selected for the FROG study. Any larger pitch amplitude would cause the wake to 
impact the horizontal tail. An example BINP8A input line is shown below. Note that 
pitch amplitude is in degrees and frequency is in rad/sec: 

&BINP8A PHIMAX=0.0, THEMAX=1.5, PSIMAX=0.0 

WTX=0.0, WTY=6.283, WTZ=0.0, &END 



Figure A.6 Pitching Motion for Obtaining the Pitch Rate Damping Derivatives. 


The number of time steps and time step interval is chosen to create a nice 
sinusoidal output through at least two cycles of motion. In the FROG study, a pitching 
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motion frequency of co=27t rad/s or one cycle/sec is utilized. Fifty time steps are chosen 
with an interval of 0.05 seconds. This combination yields 2.5 cycles of pitching motion. 
After the solution is obtained, the "total coefficient" data is extracted for plotting. The 
data can be picked out manually, or a data retrieval program can be created for the task. 
For the FROG study, the data is extracted manually and pasted into an "Excel" 
spreadsheet for plotting. In addition to the CMARC output, pitch motion phase angle and 
angle-of-attack are calculated using Equations A.33 and A.34. Lift and pitching moment 
phase angle is obtained with respect to angle-of-attack. Usually, the second motion cycle 
is used to ensure that all start up transients have subsided. Table A.5 shows 
representative FROG data for 20 time steps. 


^pitch timestep*^^ * ^pitch * ^ degrees 
a = e = sini(ftpitch)Apitch degrees 


Where: 


<j>pitch - pitching motion phase angle 

a - Angle-of-attack from the pitching motion 

#timestep - time Step 

dt - time step interval 

©pitch - pitching frequency 

Apitch - pitch amplitude 


A.33 

A.34 
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Step 

( t ) ( deg ) 

a 

CL 

CD 

CY 

C m 

C n 

c l 

0 

0 

0.00 

0.0540 

0.0273 

0.0000 

0.0828 

0.0000 

0.0000 

1 

18 

0.46 

0.4545 

0.0364 

0.0000 

- 0.0462 

0.0000 

0.0000 

2 

36 

0.88 

0.4961 

0.0344 

0.0000 

0.0303 

0.0000 

0.0000 

3 

54 

1.21 

0.5312 

0.0347 

0.0000 

0.0226 

0.0000 

0.0000 

4 

72 

1.43 

0.5491 

0.0351 

0.0000 

0.0270 

0.0000 

0.0000 

5 

90 

1.50 

0.5541 

0.0348 

0.0000 

0.0310 

0.0000 

0.0000 

6 

108 

1.43 

0.5435 

0.0338 

0.0000 

0.0447 

0.0000 

0.0000 

7 

126 

1.21 

0.5227 

0.0321 

0.0000 

0.0535 

0.0000 

0.0000 

8 

144 

0.88 

0.4925 

0.0302 

0.0000 

0.0618 

0.0000 

0.0000 

9 

162 

0.46 

0.4557 

0.0281 

0.0000 

0.0692 

0.0000 

0.0000 

10 

180 

0.00 

0.4167 

0.0262 

0.0000 

0.0727 

0.0000 

0.0000 

11 

198 

- 0.46 

0.3787 

0.0247 

0.0000 

0.0741 

0.0000 

0.0000 

12 

216 

- 0.88 

0.3457 

0.0236 

0.0000 

0.0720 

0.0000 

0.0000 

13 

234 

- 1.21 

0.3218 

0.0230 

0.0000 

0.0648 

0.0000 

0.0000 

14 

252 

- 1.43 

0.3075 

0.0227 

0.0000 

0.0585 

0.0000 

0.0000 

15 

270 

- 1.50 

0.3055 

0.0228 

0.0000 

0.0499 

0.0000 

0.0000 

16 

288 

- 1.43 

0.3158 

0.0233 

0.0000 

0.0402 

0.0000 

0.0000 

17 

306 

- 1.21 

0.3369 

0.0241 

0.0000 

0.0323 

0.0000 

0.0000 

18 

324 

- 0.88 

0.3674 

0.0254 

0.0000 

0.0246 

0.0000 

0.0000 

19 

342 

- 0.46 

0.4039 

0.0271 

0.0000 

0.0191 

0.0000 

0.0000 

20 

360 

0.00 

0.4429 

0.0290 

0.0000 

0.0153 

0.0000 

0.0000 


Table A.5 Sample Pitch Motion Data for 20 Time Steps 


The Cl and Cm responses are plotted as a function of pitch motion phase 
angle. Angle-of-attack is also plotted on the right hand vertical axis. Figure A.7 is a 
representative plot of CMARC data for the FROG UAV study. The out-of-phase 
(imaginary) portion of Cl or Cm is the combined a-dot and pitch damping contribution. 
As seen in Equations A.35 and A.36, the a-dot contribution is subtracted from the total to 
yield the pitch damping influence. The pitch rate contribution is normalized by pitch rate 
and t*. The phase angle is measured between the lift or pitching moment response and 
angle-of-attack. Equations A.35 through A.44 are used to solve for CLq and Cmq. 
Equations A.35 and A.36 assume that the a-dot and pitch rate contributions are in phase 
with each other. The a-dot contribution is small, so there isn't harm done if the two 
contributions are somewhat out-of-phase. Consequently, the a-dot contribution is 
calculated based on the maximum a-dot rate observed. 


A.35 




101 


A.36 
































J= J 


a-^q 




a+q . 


‘[CLiVcLi ‘amax*'' 


C/ 


“ * 9max * ^ 


^ J 


C 


mr 


^ ^THq * 9max * ^ 




9max “ ^max “ 


'deg 


57.3 


deg> 


0 ) 


pitch 


^rad 


* ^ rr 

' 'I^'o 


rad/s 


A.37 

A.38 

A.39 

A.40 

A.41 

A.42 

A.43 

A.44 


Where; 


[CLa-dot+q] 
[C m a-dot+q] 

I [Cl a-dot+q] 
I[C in a-dot+q] 
I[CLq] 
I[Cn,q] 

I[Cl a-dot] 

I[C 

ma- dot] 


- Amplitude of Cl response from pitching motion 

- Amplitude of Cm response from pitching motion 

- Out-of-phase Cl due to a-dot and pitch rate damping 

- Out-of-phase Cm due to a-dot and pitch rate damping 

- Cl coefficient contribution from pitch rate damping 

- Cm coefficient contribution from pitch rate damping 

- Cl coefficient contribution from a-dot damping 

- Cm coefficient contribution from a-dot damping 

- Characteristic time 


102 







NPS FROG UAV 
Oscillating Pitch Motion 
co=2pi rad/s, A=±1.5® Theta 



Figure A.7 Representative Pitch Motion Data for the FROG UAV. 


b. Sample Pitch Rate Damping Data Reduction 

Sample pitch rate damping data reduction is presented below in Equations 
A.45 through A.54. The phase angle between angle-of-attack and the Cl or Cm response 
was measured graphically by zooming the graph axes to expand the area of interest. 


20m 


t* =12.= 2 

^0 /1056 


= 0.009470 sec 
in 

sec 


9max =amax = T * pUcH mrf/s = 0.1645 rad/s 


57.3 


degy 


57.3' 


'rad 


'rad 


l\Ci^]=Ci^ *dniax *t* =1.420*0.1645*0.009470 = 0.002212 


A.45 


A.46 


A.47 
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‘[Cma ]= * “max *'" = -6264 • 0.1645 • 0.009470 = -0.009758 

1= J= * tO'234] = 0.0129 


J= sin(j»)[c„^^^ ]= sin(-l 00“) • [0.02853] = -0.02810 


'^a+g 


/[q^ J=/[C^^^^ J- /[c^^ J=0.0129 -0.002212 = 0.01069 


A.48 

A.49 

A.50 

A.51 


J- /[c„^ ]= -0.02810 - (-0.009758) = -0.01834 A.52 


Cu = 


J 


0.01069 


^ a *t 0.1645*0.009470 
v/max ‘ 


= 6.862 


A.53 


C = 


J 


-0.01834 


^ a *t 0.1645*0.009470 
ymax ‘ 


= -11.78 


A.54 
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c. 


LATERAL-DIRECTIONAL STABILITY DERIVATIVES 


This section will describe CMARC methods for the development of lateral- 
directional stability derivatives to include the static, rate damping and control power 
derivatives. The results obtained from CMARC should be spot checked against classical 
design calculations. 

Of note, CMARC uses the semi-span to normalize rolling and yawing moment 
coefficients. Most texts on stability and control, including Roskam [Ref 12] and Etkin 
[Ref 13], normalize rolling and yawing moments by span. This study will normalize roll 
and yaw moments by span. Table A.l summarizes the factors for normalizing moments 
and angular rates. All rolling and yawing moment coefficients presented in this study 
have been normalized with span by dividing the CMARC output by a factor of two. 
Table A.l also indicates the characteristic time, t*, employed for angle rate data 
reduction. 

For the lateral-directional analysis, both sides of the body must be modeled. The 
easiest way to do this is to activate the symmetric patch toggle for each patch 
(IPATSYM=1). Then, turn off symmetric calculations (RSYM=1.0). This creates 
symmetric patches around the y=0 plane, allowing CMARC to perform asymmetric 
calculations around the entire body. Processing times are significantly increased 
compared to the symmetric case. Wakes must be defined for each side of the model 
based on the new patch niunbers. In addition, don’t forget to update the adjacent patch 
numbers of the wingtips. 

The vertical tail is activated for the lateral-directional study. This requires that the 
wing and fuselage group be run separately from the horizontal and vertical stabilizer 
group. This de-coupling is required to keep the wing wake from hitting the vertical 
stabilizer. The two sets of results are summed to get the total response. Typically, the 
sidewash derivative, ds/d|3 is small, making separate solutions feasible. The main draw 
back is that separate solutions will not capture the finer interactions that a complete 
model would capture. However, many aircraft configurations will not require the de¬ 
coupling of the empennage surfaces. If the configuration permits, run the lateral- 
directional test cases as a complete airframe model. Figure A.8 shows the FROG model 
configurations used to find the lateral-directional derivatives. 
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It should be noted that the CMARC wind axes are modeled with X-aft/Z-up, vice 
X-forward/Z-down for the typical stability axes system. Positive yaw angle in CMARC 
creates positive sideslip in the stability axes system. Care must be taken to reverse the 
signs of the appropriate coefficients to convert firom the CMARC wind axes to the 
stability axis system shown in Figure A.l. 



z 

Lf . 


X 


Figure A.7 FROG UAV Lateral-Directional Models. 

1. Static Lateral Directional Stability Derivatives 

a. Static Lateral-Directional Stability Equations 

Three lateral-directional stability derivatives are measured with the first 
two runs of the CMARC model. The model is first analyzed at the estimated trim 
condition with zero sideslip. The baseline solution should be compared to the 
longitudinal baseline. If the lift coefficient differs significantly, the wake modification 
assumptions may need to be revisited. The model is then checked for lateral-directional 
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balance at zero sideslip angle (yaw angle =0). The side force, rolling and yawing 
coefficients should be zero when a trial run is performed at zero sideslip. If lateral- 
directional forces or moments are present, the model and wake geometry should be 
checked for symmetry. Visualization with LOFTSMAN will assist in spotting the 
problem. 

For the FROG case, art„=0° and P=0° is selected for the baseline. A 
second CMARC run is then conducted with yaw angle incremented by one or two 
degrees. Cy, C, and C„ are then extracted manually fi-om the data files. The slope of Cy, 
C, and C„ versus sideslip provides the Cyp, C,p and C„p longitudinal derivatives. 
Additional solutions can be obtained to check for slope linearity. Equations A.55 through 
A.57 are used for these calculations: 

^ Cl . 180 ,. 

Cia = —— * -per radian A.56 

P A/3° 7t 

C 180 

=——* -per radian A.57 

A/3° n 

b. Sample Static Lateral-Directional Data Reduction 

Table A.6 presents CMARC solutions for the FROG UAV at two angles- 
of-attack. Data is for the summed contribution of the wing/fuselage and tail group 
models. Roll and yaw moments are already divided by a factor of two. This normalizes 
by span to compensate for the CMARC output, which is normalized by semi-span. 
Sample calculations for obtaining the static derivatives are demonstrated below. The 
calculations are easily implemented in a spreadsheet or with a MATLAB script. 


RUN# 

3 

CL 

CD 

CY 

C m 

C n 

C 1 

1 

0° 

0.4432 

0.0145 

0.0000 

0.0048 

0.0000 

0.0000 

2 

2° 

0.4452 

0.0148 

-0.0087 

0.0038 

0.0022 

-0.0022 


Table A.6 FROG UAV Static Lateral-Directional Stability Data at atrio,=0°. 
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Sample Calculations: 


^ Cy *180 -0.0087,180 deg .. 

Cvr, = —^ *-=--*- — = -0.2493 per radian 


'Y/S 




TT 


n rad 



Cl * 180 _ -0.0022 * 180 deg 
Ayff° n 2 7Z rad 


= 0.0630 per radian 


A.55 


A.56 


Ay9 


O 


n 


0.0022 , 180 deg 
2 ^ rad 


= 0.0630 per radian 


A.57 


3. Lateral-Directional Control Power Stability Derivatives 

a. Lateral-Directional Control Power Equations 

The rudder control power derivatives are obtained by substituting a 0° 
deflection vertical stabilizer patch for one with positive rudder deflection. Only one run 
is required. For the FROG UAV study, +5° (TEL) deflection is used. The difference 
between the trim condition and the deflected value is divided by the rudder deflection as 
shown below; 




C 





^Y5a2 

-^Ysai 

o 

00 

_* 


-Sai 

n 



*180 

5a2 - 

-5a\ 

n 

l^«&2 

-C 

1.180 

&2 

- 5ai 

n 

■ 

■ 1 


Sr2 - 

-Srj 

n 


per 


per 


per 


per 


rad 


rad 


rad 


rad 


A.58 


A.59 


A.60 


A.61 
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A.62 


^^dri ^^5r\ i 180 , 

=- - - — *- per rad 

Sr2 - 5ri 

180 j 

=-=- — *- per rad 

Sr2 - Sr^ 


A.63 


b. Sample Lateral-Directional Control Power Data Reduction 


Tables A.7 and A.8 present CMARC solutions for aileron and rudder 
deflections. Sample calculations for obtaining aileron and rudder control power are 
demonstrated below. The calculations are easily implemented in a spreadsheet or with a 
MATLAB script. 


'YSa 


^ ^YSa2 180 _ (-0.0018-0) ^ 180 _ Q 

Sa2 - <5hi 7C (5 - O) K 


0103 


A.64 


Q, 


Sa 


^ ~^/&ii ^ l80 ^ (0.01695-0) ^ 180 _Q 


6a2 -da^ 


n 


(5-0) 


1943 


K 


^ 180 , (-000105-0) ^ 180 _ ^ 

5a2 — Soi n (^ ~ O) n 


A.65 


A.66 


'YSr 


-Cf&,) .180 (0.0081-0) .180 


5r2 - <5irj 


n 


(5-0) 


n 


A.67 


Q, 


8r 


_ ^ hr 2 ~^lSri 180 (0.00035-0) ^ 180 


Sr2 - Sri 


n 


( 5 - 0 ) 


0040 


n 


A.68 


6r2-6ri 


nsri i^l80 _ (-0.00395-0)^180 _ 


TT 


( 5 - 0 ) 


= -0.0453 


7t 


A.69 
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RUN# 

6a 

CL 

CD 

CY 

C m 

C n 

C 1 

1 

0° 

0.4259 

0.0170 

0.0000 

0.0104 

0.0000 

0.0000 

2 

5° 

0.4005 

0.0173 

-0.0018 

0.0199 

-0.00105 

0.01695 


Table A.7 FROG UAV Aileron Control Power Data at =0°. 


RUN# 

5r 

CL 

CD 

CY 

C m 

C n 

C 1 

1 

0° 

0.4259 

0.0170 

0.0000 

0.0104 

0.0000 

0.0000 

2 

5° 

0.4258 

0.0170 

0.0081 

0.0107 

-0.00395 

0.00035 


Table A.8 FROG UAV Rudder Control Power Data at a^„=0°. 

3. Yaw Rate Damping Derivatives 

a. Yaw Rate Derivative Methods and Equations 

Only one motion is required for the yaw rate terms. The |3-dot terms are 
generally considered negligible. Therefore, a sideways plunging motion is not required. 
The yaw rate terms are yielded directly from an oscillating yawing motion as depicted in 
Figure A.8. Yawing motion data for the FROG is gathered at a frequency of 1% rad/s, 
which equates to a reduced frequency of k=0.369 for this configuration and trim airspeed. 
Oscillating yaw motion is controlled with the CMARC BINP8A input file line. An 
amplitude of ±2° is selected for the FROG study. An example BINP8A input line is 
shown below. Note that yaw amplitude is in degrees and frequency is in rad/sec: 

&BINP8A PHIMAX=0.0, THEMAX=1.5, PSIMAX=0.0 

WTX=0.0, WTY=6.283, WTZ=0.0, &END 


no 
































Figure A.8 Yawing Motion for Obtaining the Yaw Rate Derivatives. 


The number of time steps and time step interval is chosen to create a nice 
sinusoidal output through at least two cycles of motion. In the FROG study, a yawing 
motion frequency of co=2n rad/s or one cycle/sec is utilized. Fifty time steps are chosen 
with an interval of 0.05 seconds, which creates 2.5 cycles of yawing motion. After the 
solution is obtained, the "total coefficient" data is extracted for plotting. The data can be 
picked out manually, or a data retrieval program can be created for the task. For the 
FROG study, the data is extracted manually and pasted into an "Excel" spreadsheet for 
plotting. As discussed earlier, the FROG model was split into wing/fuselage and tail 
groups for the lateral-directional study. Both sets of data are pasted into the Excel 
spreadsheet and summed. The total value is then plotted. In addition, yawing motion 
phase and sideslip angles are also calculated using Equations A.70 and A.71. Yawing 
motion phase angle is plotted on the x-axis and sideslip is plotted on the right hand axis. 
Table A.9 shows representative FROG yaw rate data for 20 time steps. 


fjy -u . *rrT * 

ryaw ^timestep ^ yaw 


180 


K 


degrees 


A.70 
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P = * Ayaw degrees 


A.71 


Where: 


*))yaw - yawing motion phase angle 

P - Sideslip angle from yawing motion 

#dmes.ep ' time Step 

dt - time step interval 

- Yawing frequency 

- Yaw angle amplitude 


Step 

(|) ( deg ) 

p 

CL 

CD 

CY 

C m 

C n 

c l 

0 

0 

0.00 

0.0534 

0.0111 

- 0.0005 

0.085 

0.0003 

- 0.0002 

1 

18 

0.62 

0.4004 

0.0195 

- 0.007 

- 0.0107 

0.0023 

- 0.0023 

2 

36 

1.18 

0.4329 

0.0161 

- 0.0087 

- 0.0009 

0.0026 

- 0.0026 

3 

54 

1.62 

0.4388 

0.0152 

- 0.0098 

0.0019 

0.0027 

- 0.0027 

4 

72 

1.90 

0.4416 

0.0149 

- 0.0098 

0.0027 

0.0026 

- 0.0026 

5 

90 

2.00 

0.4428 

0.0149 

- 0.0089 

0.0033 

0.0022 

- 0.0022 

6 

108 

1.90 

0.4435 

0.0147 

- 0.0071 

0.0039 

0.0016 

- 0.0016 

7 

126 

1.62 

0,444 

0.0146 

- 0.0046 

0.0044 

0.0008 

- 0.0008 

8 

144 

1.18 

0.4444 

0.0146 

- 0.0017 

0.0049 

0.0000 

0.0000 

9 

162 

0.62 

0.4447 

0.0144 

0.0014 

0.0051 

- 0.0009 

0.0009 

10 

180 

0.00 

0.445 

0.0144 

0.0043 

0.0048 

- 0.0016 

0.0016 

11 

198 

- 0.62 

0.4452 

0.0145 

0.0069 

0.0044 

- 0.0022 

0.0022 

12 

216 

- 1.18 

0,4453 

0.0146 

0.0087 

0.004 

- 0.0026 

0.0026 

13 

234 

- 1.62 

0.4452 

0.0147 

0.0098 

0.0035 

- 0.0027 

0.0027 

14 

252 

- 1.90 

0.4451 

0.0147 

0.0098 

0.0034 

- 0.0026 

0.0025 

15 

270 

- 2.00 

0.4449 

0.0148 

0.0089 

0.0036 

- 0.0022 

0.0021 

16 

288 

- 1.90 

0.4448 

0.0146 

0.0071 

0.0041 

- 0.0016 

0.0015 

17 

306 

- 1.62 

0.4449 

0.0146 

0.0046 

0.0045 

- 0.0008 

0.0008 

18 

324 

- 1.17 

0.445 

0.0145 

0.0017 

0.005 

0.0000 

- 0.0001 

19 

342 

- 0.62 

0.4452 

0.0144 

- 0.0014 

0.0051 

0.0009 

- 0.0009 

20 

360 

0.00 

0.4454 

0.0144 

- 0.0044 

0.0049 

0.0016 

- 0.0016 


Table A.9 Sample FROG UAV Yawing Motion Data for 20 Time Steps 


The Cy, C, and responses are plotted as a function of yawing phase 
angle. Sideslip is also plotted on the right hand vertical axis. Figure A.9 is a 
representative plot of CMARC data for the FROG UAV study. The out-of-phase 
(imaginary) portion of Cy, C, or is due to the yaw rate damping contribution. The p- 
dot contribution is considered small and is ignored. The phase angle is measitred 
between the coefficient response and sideslip angle. In this case, the zero crossing 


112 

































method of measuring phase angle proves to the easiest and most accurate. The phase 
angle is simply measured from the parameter to sideslip angle where they cross the x-axis 
in the same direction. The yaw rate damping contribution is normalized by maximum 
yaw rate and t*. Equations A.72 through A.76 are used to solve for Cy,, C,, and C„. 



I[Cy\ 


f ^ t 

'^max * 


* ’ 


where 


I[Cy] = sini^ya.^ \Cy ] 


Cir " ^[Q ] = sin((z5^aw \Cl ] 

Q, yyhere l[Cn]-sm{^ya^\Cn\ 

^max 


Vdeg 


'max 


57.3 


deg/ 


0) 


yaw 


'rad 


* /2 
t =-^ sec 

Ur, 


rad/s 


A.72 

A.73 

A.74 

A.75 

A.76 


Where: 

[Cy] - Amplitude of Cy response from yawing motion 

[C,] - Amplitude of C] response from yawing motion 

[CJ - Amplitude of C„ response from yawing motion 

I[Cy] - Cy coefficient contribution from yaw rate damping 

I[C|] - C, coefficient contribution from yaw rate damping 

I[C„] - C„ coefficient contribution from yaw rate damping 

t* - Characteristic time 
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Figure A.9 Representative Yawing Motion Plot for the FROG UAV. 


b. Sample Yaw Rate Damping Data Reduction 

Sample data reduction is presented below in Equations A.28 through A.32. 
The phase angle between the response and the sideslip angle is measured graphically by 
zooming in on the x-axis zero crossing. 


b/ \2Ain/ 

t* = -^ = 0.0587 sec A.77 

Uq 1056 ^% 



CO 


yaw 



'max ” 


0.2193 rad/s 


A.78 

































































A.79 


Q ACy] _ sin(^jvawjCy] _ sin(l54°X0.0099] 


'max 




r ^ t 
'max ' 


0.2193*0.0587 


= 0.3371 


Q /[Q ] J[C/ ] sin(l 44°X0.00265] 


V ^ t V ^ t 

'max ' 'max ^ 


0.2193*0.0587 


= 0.1210 


C„ 


ACn] _ _ sin(-36°X0.00265] 

* * ____ ^ 




'max 


*t 


0.2193*0.0587 


1210 


A.80 


A.81 


4. Roll Rate Damping Derivatives 

CL Roll Rate Derivative Methods and Equations 

Gathering roll damping data is straightforward. Unlike the pitch or yaw 
rate terms, there is no change in angle-of-attack or sideslip with the rolling motion. 
Therefore, roll damping effects can be measured with pure rolling motion around the x- 
axis. As with the yaw damping case, the wing/fuselage and tail models are run 
separately. The results are then summed. A rigid wake seems to work well, which 
provides a solution after just a few time steps. 

Roll damping data for the FROG is gathered at a 20 deg/sec roll rate. 
Initially, a 5 deg/sec roll rate was selected, but the yawing moment coefficient was low 
enough that numerical resolution would effect the results. The higher rate provided a 
sufficiently large response for all terms. Pure rolling motion is controlled by the 
PHIDOT term in the CMARC BINP8 input file line. An example BINP8 input line is 
shown below. Note that roll rate is in deg/sec: 

&BINP8 ALDEG=0.0, YAWDEG=0.0, PHIDOT=20.0 

THEDOT=0.0, PSIDOT=0.0, &END 

The roll rate damping terms are obtained by normalizing the CMARC 
output by roll rate, p, and characteristic time, t*. Equations A.82 through A.85 are used 
for these calculations. 
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A.82 




Q 


p*t 


* 


C =-^ 


p*t 


* 

t 



sec 


Where; 

[Cy] - Amplitude of Cy response from roll rate 

[C,] - Amplitude of C, response from roll rate 

[C„] - Amplitude of C„ response from roll rate 

p - Roll rate 

t* - Characteristic time 


A.83 


A.84 


A.85 
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d. Sample Roll Rate Damping Data Reduction 


Tables A. 10 presents the CMARC roll damping solution for the FROG 
UAV. The data is the summed contributions from the wing/fliselage and tail surface 
models. Sample roll damping calculations presented below in Equations A. 86 through 
A.89. The calculations are easily implemented in a spreadsheet or with a MATLAB 
script. 


t* =^= . = 0.0587 sec 

Uq \Q56W 


A.86 


Cy^ = 


Cu = 





0.0010 

p*t 

20 


^ 0.0587 sec 



sec 

180 rad 

Cl 



-0.00925 

p*t* 

20 

deg* 

^ 0.0587sec 
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-0.00045 

p*t* 

20 

deg* 

^ *0.0587 sec 



sec 

180 rad 


= 0.0488 
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A.89 


RUN# 

P 
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C m 
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1 

0 deg/s 
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20 deg/s 

0.4435 
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0.0047 

-0.00045 
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Table A.IO FROG UAV Roll Damping Data at and P=0®. 
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&BPNODE TN0DE=3, TNPC=0, T1NTC=0. «.END 4BPNODE TNODE=3. TNPC=0, T1NTC=0. iENH 

&SECT1 STX = 0.0, STY = 0.0. STZ = 0.0. SCALE=1.0, ALF=0.0, THETA=0.0, INM0DE=4. 4SECT1 STX = 0.0. STY = 0.0. STZ=0.0. SCALE«=1.0. ALF=0.0, THETA^O.O. INMODE 

TNODSaO. TNPS = 0. T1NTS = 0, iEND TNODScO, TNPS=0, T1NTS«=0, &END 

44.4894 53.5359 11.5387 43.8419 57.0000 11.5897 
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APPENDIX D 


MATLAB SCRIPT FOR DYNAMIC MODE ANALYSIS 
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% froguav.m: This MATLAB script calculates FROG UAV dynamics. Eigenvalue modal analysis 
% and response to control input is obtained for three aerodynamic models. 

% 

clear 


% Calculate trim lift and drag coefficients 

% FROG Trim condition 

W = 67.73; 
m = W/32.2; 

S = 2530/144; 

V = 88; 

M = V/1118; 
cbar = 20/12; 
b = 126.5/12; 
rho = 0.002327; 
cruise_HP = 5; 
etaP = 0.35; 

L_D=7; 

CLO = 0.4295; 

Ixx = 12.52; lyy = 
q = 0.5*rho*(V)^2; 

CL = W/(q*S); 

% Input the non-dimensional stab derivatives: 

% ★****★**★*★★★*★★★*★★*★*★***★******★*★★***★★**★**★★*★**★★*************★★********★******** 
for i = 1:3 


% lbs 
% slugs 

% ft"2, 2530 in^2 
% ft/s = 1056 in/s 
% trim mach number 
% ft = 20 in 
% ft = 124 in 

% slugs/ft^3 800 ft MSL std day 


8.43; Izz = 18.55; lxz=0; %slugs-ft'^2 


if i==l 


source=' CMTUIC ’ 

% Longitudinal Stability Derivatives 

D_L_D=W/L_D;CD_L_D = D_L_D/(q*S) ; % Drag from Lift to Drag Ratio 

Tr=cruise_HP*etaP*550/V; CD_Tr=Tr/(q*S); % Drag from thrust required calculation 

CD=(CD_L_D +CD_Tr)/2; % CD averaged from both methods 

CLalpha = 4.845; CDalpha = 0.2664; CMalpha = -0.4126; 

CLalphadot = 1.420; CMalphadot = -6.264; CDalphadot =0; % Assumed to be small 

alphatrim = 0; % radians 

thetaO = alphatrim; 

CLmach = 0; CDmach = 0; CMmach = 0; 

CLq = 6.82; CMq = -11.78; CDq =0; % Assumed to be small 

CLde = 0.43 78; CMde = -1.199; CDde = 0.0092; % Induced drag contribution only 

% Lateral - Directional Stability Derivatives 

CYbeta = -0.2493; Clbeta = -0.0630;Cnbeta = 0.0630; 

Clp = -0.4514; Cnp = -0.0220; CYp = 0.0488; 

Clr = 0.1210; Cnr = -0.1210; CYr = 0.3370; 

Cldr = 0.0040; Cndr = -0.0453; CYdr = 0.0928; 

Clda = 0.1943; Cnda = -0.0121; CYda = -0.0206; 

elseif i==2 

source='Classical_Analysis' 

% Longitudinal Stability Derivatives 
D_L_D=W/L_D; CD_L_D = D_L_D/(q*S); 

Tr=cruise_HP*etaP*550/V; CD_Tr=Tr/(q*S); 

CD=(CD_L_D +CD_Tr)/2; 

CLalpha = 4.82; CDalpha = 0.253; CMalpha = -0.70 % d_epsilon/d_alpha = 0.40 

CLalphadot = 1.56; CMalphadot = -4.14; CDalphadot =0; % Assumed to be small 

alphatrim = 0; % radians 

thetaO = alphatrim; 

CLmach = 0; CDmach = 0; CMmach = 0; 

CLq = 3.89; CMq = -11.39; CDq =0; % Assumed to be small 

CLde = 0.39; CMde = -1.04; CDde = 0.00; 


% Drag from Lift to Drag Ratio 
% Drag from thrust required calculation 
% CD averaged from both methods 
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% Lateral - Directional Stability Derivatives 
CYbeta = -0.511; Clbeta = -0.055;Cnbeta = 0.051; 

Clp = -0.300; Cnp = -0.072; CYp =0; % Assumed to be small 

Clr = 0.168; Cnr = -0.0762; CYr = 0.140; 

Cldr = 0.0056; Cndr = -0.0341; CYdr = 0.081; 

Clda = 0.213; Cnda = -0.0236; CYda = -0.00; 

else 


source='Parameter_Analysis' 

% Longitudinal Stability Derivatives 
CL = 0-4295; CD = 0.0614; 

CLalpha = 4.0907; CDalpha = 0.23; CMalpha = -0.4174; 

CLalphadot = 1.3877; CMalphadot = -3.7115; CDalphadot =0; % Assumed to be small 

alphatrim =0; % radians 

thetaO = alphatrim; 

CLmach = 0; CDmach = 0; CMmach = 0; 

CLq = 3.35; CMq = -8.8818; CDq =0; % Assumed to be small 

CLde = 1.1249; CMde = -1.6208; CDde = 0.0676; 

% Lateral - Directional Stability Derivatives 

CYbeta = -0.9867; Clbeta = -0.0942;Cnbeta = 0.1755; 

Clp = -0.4483; Cnp = -0.1077; CYp =0; % Assumed to be small 

Clr = 0.2078; Cnr = -0.1212; CYr = 0.1096; 

Cldr = 0.0004; Cndr = -0.0785; CYdr = 0.0926; 

Clda ^ 0.2387; Cnda = -0.0261; CYda = 0; 


end 


% Calculate the Dimensional Stability Derivatives 


Xu = -q*S/(m*V) * (2*CD+M*CDmach) ; % Xu - 1/sec , assumes dCD/dM=0 (ie no Mach effect) 

Xalpha = q*S/m*(CL-CDalpha); % Xalpha - ft/sec^2 

Xalphadot = -q*S/ra*cbar/(2*V)*CDalphadot; % Xalphadot - ft/sec 

Xq = -q*S/m*cbar/(2*V)*CDq; % Xq - ft/sec 

Xde = -q*S/m*CDde; % Xde - ft/sec^2 


Zu = -q*S/(m*V)*(2*CL+M*CLmach); 
Zalpha = -q*S/m*(CD+CLalpha) 


% Zu - 1/sec , assumes dCL/dM=0 (ie no Mach effect) 
% Zalpha - ft/sec^2 


Zalphadot = -q*S/m*cbar/(2*V)*CLalphadot; 
Zq = -q*S/m*cbar/(2*V)*CLq ; 

Zde = -q*S/m*CLde; 


% Zalphadot - ft/sec 
% Zq - ft/sec 
% Zde - ft/sec^2 


Mu = q*S*cbar/(Iyy*V)*CMmach; 

Malpha = q*S*cbar/Iyy*CMalpha; 

Malphadot = q*S*cbar/Iyy*cbar/(2*V)*CMalphadot ; 
Mq = q*S*cbar/Iyy*cbar/(2*V)*CMq; 

Mde = q*S*cbar/Iyy*CMde; 


% Mu - ft/sec 
% Malpha - l/sec"2 
% Malphadot - 1/sec 
% Mq - 1/sec 
% Mde - l/sec^2 


% Linearized Longitudinal 4x4 Plant Matrix 
% Form the An Bn and In plant matrices 

An = [V*Xu Xalpha 0 -32.17*cos(thetaO); V*Zu Zalpha (V+Zq) -32.17*sin(thetaO); 
V*Mu Malpha Mq 0; 0010); 

Bn = [Xde Zde Mde 0)'; 

In = [V 0 0 0; 0 (V-Zalphadot) 0 0; 0 -Malphadot 1 0; 0 0 0 1] ; 

A = inv{In)*An; 

B = inv(In)*Bn; 

% Find Short and Long Period Natural Frequency and damping 
out=source 

out=’Longitudinal mode E-values' 

P = poly(A); 
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R = roots(P) 
out='Short Period' 
phi = angle(R); 

Zl(i) = (l/(H-tan(phi(l) ) "2) ) ".5 
Wnl(i) = real(R(l))/(-Z1 (i)) 

Wdl(i) = (l-Zl(i)"2)".5*Wnl(i); 

Tdl(i) = 2*pi/Wdl(i); 

out=’Long Period' 

Z3(i) = (1/(l+tan(phi(3))"2))".5 
Wn3(i) = real(R(3))/(-Z3(i)) 

Wd3(i) = (l-Z3(i)"2)".5*Wn3(i); 

Td3(i) = 2*pi/Wd3(i); 

out='Step response' 

Cl=[l 0 0 0;0 1 0 0/0 0 1 0; 0 0 0 1]; D1=0; 
U=zeros{1001,1) ; 

U(101:1001)=-2/57.3; 
t=(0:.01:10); 

SYS=SS(A,B,C1,D1); 

[Y,T]=lsim(SYS,U,t); 
tl=t; 

de_s=U*57.3; 

alpha_s(:,i)=Y(:,2)*57.3; 
pr_s(:,i)=Y(:,3)*57.3; 
theta_s(: , i)=Y(:,4)*57.3; 

Cl=[l 0 0 0;0 1 0 0;0 010; 0001]; D1=0; 
U=zeros(1001,1) ; 

U(101:160)=-5/57.3; 

U(161:220)=5/57.3; 
t=(0:.01:10); 

SYS=Ss(A,B,Cl,Dl); 

[Y,T]=lsim(SYS,U,t); 
tl=t; 

de_d=U*57.3; 

alpha_d(: ,i)=Y(:,2)*57.3 ; 
pr_d ( : , i) =Y ( : , 3) * 57 . 3 ; 
theta_d(:,i)=Y(:,4)*57.3; 


% Lateral-Directional 4x4 Plant Matrix 


% Ybeta - ft/sec^2 
% Yp - ft/sec 
% Yr - ft/sec 
% Yda - ft/sec^2 
% Ydr - ft/sec^2 

% Lbeta - l/sec^2 
% Lp - 1/sec 
% Lr - 1/sec 
% Lda - l/sec^2 
% Ldr - l/sec^2 

% Nbeta - l/sec^2 
% Np - 1/sec 
% Nr - 1/sec 
% Nda - l/sec^2 
% Ndr - l/sec^2 


% Linearized Longitudinal 4X4 Plant Matrix 
source 


% Calculate the Dimensional Stability Derivatives 

Ybeta = q*S/m*(CYbeta); 

Yp = q*S/m*b/(2*V)*CYp; 

Yr = q*S/m*b/(2*V)*CYr; 

Yda = q*S/m*CYda; 

Ydr = q*S/m*CYdr; 

Lbeta = q*S*b/Ixx*(Clbeta); 

Lp = q*S*b/Ixx*b/(2*V)*Clp; 

Lr = q*S*b/Ixx*b/(2*V)*Clr; 

Lda = q*S*b/Ixx*Clda; 

Ldr = q*S*b/Ixx*Cldr; 

Nbeta = q*S*b/lzz*(Cnbeta); 

Np = q*S*b/Izz*b/(2*V)*Cnp; 

Nr = q*S*b/Izz*b/(2*V)*Cnr; 

Nda = q*S*b/lzz*Cnda; 

Ndr = q*S*b/Izz*Cndr; 
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longTnodes=' Lateral Directional 4x4 Plant Matrix' 

% Form the An Bn and In plant matrices 

Cn = [Ybeta Yp 32.17*cos (thetaO) (Yr-V) ; Lbeta Lp 0 Lr; 0 10 0; Nbeta Np 0 Nr] 
Dn = [Yda Lda 0 Nda; Ydr Ldr 0 Ndr]'; 

Iln = [V 0 0 0; 010 -Ixz/Ixx; 0 0 1 0; 0 -Ixz/Izz 0 1]; 

C = inv(IIn)*Cn; 

D = inv(Iln)*Dn; 

% Find Dutch roll, roll and spiral Natural Frequency and damping 
out='Lateral Directional mode E-values' 

PI = poly(C); 

R1 = roots(PI) 
if i==l 

out = ' Dut ch Rol 1 Mode ' 
phi = angle(Rl); 

Z2(i) = (l/(l+tan(phi(2))"2))".5 
Wn2(i) = real(R1(2))/(-22(i)) 

Wd2(i) = (1-22(i)"2)".5*Wn2(i); 

Td2(i) = 2*pi/Wd2(i); 

out=’Roll Mode' 
roll(i) = (RKD) 
out='Spiral Mode' 

Spiral(i) = (Rl(4)); 

elseif i==2 

out='Dutch Roll Mode' 

phi = angle(Rl); 

22(i) = (l/(l+tan(phi(l))"2))".5 
Wn2(i) = real(Rl(1))/(-22(i)) 

Wd2(i) = (l-Z2(i)^2)".5*Wn2(i); 

Td2(i) = 2*pi/Wd2(i); 
out='Roll Mode' 
roll(i) = (Rl(3)) 
out='Spiral Mode' 

Spiral(i) = (Rl(4)); 

else 

out='Dutch Roll Mode' 
phi = angle(Rl); 

Z2(i) = (l/(l+tan(phi(l))^2))".5 
Wn2(i) = real(Rl(1))/(-22(i)) 

Wd2(i) = (1-22(i)"2)".5*Wn2(i); 

Td2(i) = 2*pi/Wd2(i); 
out='Roll Mode' 
roll(i) = (Rl(3)) 
out=’Spiral Mode' 

Spiral(i) = (Rl(4)) 


end 

% Aileron step 

Cl=[l 0 0 0;0 1 0 0;0 0 1 0; 0 0 0 1]; D1=0; 
U=zeros(1001,2); 

U(101:1001,l)=2/57.3; 
t=(0: .01:10) ; 

SYS=ss(C,D,Cl,Dl); 

[Y,T]=lsim(SYS,U,t); 
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tld=t; 

da_sa=U*57.3; 
rr_sa(:,i) =Y{:,2)*57.3; 
phi_sa(: , i)=Y{:,3)*57.3; 
beta_sa(:,i)=Y(:,l)*57.3; 

% Aileron doublet 


Cl=[l 0 0 0;0 1 0 0;0 0 1 0; 0 0 0 1]; D1=0; 
U=zeros{1001,2); 

U(101:175,l)=-5/57.3; 

U{176:250,1)=5/57.3; 
t=(0:.01:10); 

SYS=SS(C,D,C1,D1); 

[Y,T]=lsim(SYS,U,t); 
tld=t; 

da_da=U*57.3; 
rr_da(:,i)=Y(:,2)*57.3; 
phi_da(:,i)=Y(:,3)*57.3; 
beta_da{:,i)=Y(:,1)*57,3; 


% Rudder step 

Cl=[l 0 0 0;0 1 0 0;0 0 1 0; 0 0 0 1]; D1=0; 
U=zeros(101,2); 

U (101 : 1001 , 2 ) =2/57 . 3 ; 
t=(0:.01:10); 

SYS=ss(C,D,Cl,Dl); 

[Y,T]=lsitn(SYS,U,t) ; 
tld=t; 

dr_sr=U*57.3; 
beta_sr(:,i)=Y(:,1)*57.3; 
lr_sr(;,i)=Y(:,2)*57.3; 
rr_sr(:,i) =Y(: , 4)*57.3 ; 

% Rudder doublet 

Cl=[l 0 0 0;0 1 0 0;0 010; 0001]; D1=0; 
U=zeros(1001 , 2 ); 

U(10l:175,2)=-5/57.3; 

U(176:250,2)=5/57.3; 
t=(0: .01:10) ; 

SYS=ss(C,D,Cl,Dl); 

[Y,T]=lsim(SYS,U,t) ; 
tld=t; 

dr_dr=U*57.3; 
beta_dr{:,i)=Y(:,1)*57.3; 
lr_dr(:,i)=Y(:,2)*57.3; 
rr_dr(:,i)=Y(:,4)*57.3; 


end 

%****★★***★**★★★★★*★*★★★★★★★★★★★★**★**★*★*★**★★★★★****★*★★******★★*★***★★★*★*****★★★★★★★*★********★ 
% Plot Responses 
% Plot Elevator Step Response 
figure 

subplot(4,1,1), plot(tl,de_s),title('LONGITUDINAL RESPONSE TO A 2 DEGREE ELEVATOR STEP 
INPUT'),ylabel(’Elevator (deg)'),axis([0 5 -4 4]),grid 
subplot(4,1,2) , plot(tl,alpha_s (:,1),'-,',tl,alpha_s(:,2), '--',tl,alpha_s(:,3),'-') ,ylabel('A.O.A 
(deg)'),axis{[0 5 -5 10]),grid 

subplot(4,1,3), plot(tl,pr_s (:,1), '-.',tl,pr_s{:,2,tl,pr_s{:,3), '-'),ylabel ( 'q 
(deg/s) '),axis( [0 5 -20 20]),grid 

subplot(4,1,4), plot(tl,theta_s (:,1)tl,theta_s(:,2)tl,theta_s(:,3), '-'),ylabel('Theta 
(deg)'),grid 

xlabel('Time (sec)') ,axis([0 5 -10 50]), pause, 
legend('CMARC','Classical','Parameter Est ',0) 

% Plot Elevator Doublet Response 
figure 
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subplot(4,1,1) , plot(tl,de_d) ,title('LONGITUDINAL RESPONSE TO A 5 DEGREE ELEVATOR 
DOUBLET'),ylabel('Elevator (deg)'),axis([0 5 -10 10)),grid 
subplot(4,1,2) , plot(tl,alpha_dtl,alpha_d (:,2) , ',tl,alpha_d(:,3), '-'),ylabel('A.O.A 
(deg)'),axis([0 5 -12 12]),grid 

subplot(4,1,3), plot(tl,pr_d(tl,pr_d(:,2),',tl,pr_d(:,3),'-'),ylabel('q 
(deg/s)'),axis([0 5 -60 50)),grid 

subplot (4,1,4), plot (tl, theta__d (:, 1) , ' - . ' , tl, theta_d ( : , 2) , ' - - ' , tl, theta_d( : , 3) , ' - ') ,ylabel ( 'Theta 
(deg)'),grid 

xlabel('Time (sec)'),,axis([0 5 -10 30)),pause 
legend('CMARC,'Classical','Parameter Est ',0) 

% Plot Aileron Step Response 
figure 

subplot (4,1,1) , plot (tld, da__sa (:, 1) ), title ('LATERAL-DIRECTIONAL RESPONSE TO A 2 DEGREE AILERON 
STEP INPUT'),ylabel('Aileron (deg) '),axis( [0 5 -4 4)),grid 
subplot(4,1,2), plot(tld,rr_sa(tld,rr_sa(:,2) , ',tld,rr_sa(:,3), '-'),ylabel('p 
(deg/s)'),axis([0 5 -10 25]),grid 

subplot(4,1,3) , plot(tld,phi_sa(:,1), '-.',tld,phi_sa(:,2), '',tld,phi_sa(:,3), '-'),ylabel('Bank 
Angle (deg)'),axis([0 5 -10 70]),grid 

subplot(4,1,4) , plot(tld,beta_sa(:,1), '-.',tld,beta_sa(;,2),',tld,beta_sa(:,3), '- 
'),ylabel('Beta (deg)'),grid 

xlabel('Time (sec)') ,axis([0 5 -5 10)), pause, 
legend('CMARC ,'Classical','Parameter Est ',0) 

% Plot Aileron Doublet Response 
figure 

subplot(4,1,1) , plot(tld,da_da(:,1)),title('LATERAL-DIRECTIONAL RESPONSE TO A 5 DEGREE AILERON 
DOUBLET'),ylabel('Aileron (deg) '),axis( [0 8 -10 10]),grid 
subplot(4,1,2) , plot(tld,rr_da(:,1)tld,rr_da(;,2), ',tld,rr_da(;,3), '-'),ylabel('p 
(deg/s)'),axis([0 8 -50 60)),grid 

subplot(4,1,3) , plot(tld,phi_da(:,!),'-.',tld,phi_da (:,2), ',tld,phi_da(:,3),'-'),ylabel('Bank 
Angle (deg) '),axis( [0 8 -30 10)),grid 

subplot (4,1,4), plot (tld, beta_da (:,!),'-.', tld,beta__da ( : , 2) , ' - - ' , tld, beta__da (: , 3) , ' - 
'),ylabel('Beta (deg)'),grid 

xlabel('Time (sec)'),axis([0 8 -10 15)), pause, 
legend('CMARC ,'Classical','Parameter Est ',0) 

% Plot Rudder Step Response 
figure 

subplot (4,1,1) , plot (tld,dr__sr( : ,2) ), title ('LATERAL-DIRECTIONAL RESPONSE TO A 2 DEGREE RUDDER STEP 
INPUT'),ylabel('Rudder (deg) '),axis([0 8 -4 4)),grid 
subplot(4,1,2), plot(tld,beta_sr(tld,beta_sr(:,2),'',tld,beta_sr(:,3),'- 
'),ylabel('Beta (deg)'),axis([0 8 -1 3]),grid 

subplot (4,1,3) , plot (tld, rr_sr ( : , 1) , ' - . ' , tld, rr__sr ( ; , 2) , ' , tld, rr_sr ( : , 3) , ' - ' ) , ylabel ( ' r 

(deg/s)'),axis([0 8 -10 5]),grid 

subplot(4,1,4), plot(tld,lr_sr(:,1)tld,lr_sr(:,2),'--',tld,lr_sr(:,3),'-'),ylabel('p 
(deg/s)'),grid 

xlabel('Time (sec)'),axis([0 8 -5 10)), pause, 
legend('CMARC ,'Classical','Parameter Est ',0) 

% Plot Rudder Doublet Response 
figure 

subplot (4,1,1) , plot (tld, dr__dr(:, 2) ), title ('LATERAL-DIRECTIONAL RESPONSE TO A 5 DEGREE RUDDER 
DOUBLET'),ylabel('Rudder (deg) '),axis( [0 8 -10 10)),grid 
subplot(4,1,2) , plot(tld,beta_dr(tld,beta_dr(:,2), ',tld,beta_dr(:,3), '- 
'),ylabel('Beta (deg)'),axis([0 8 -10 10)),grid 

subplot (4,1,3) , plot (tld, rr_dr (:,!),'-.', tld, rr__dr ( : , 2) , ' - - ' , tld, rr_dr ( : , 3) , ' - ' ) ,ylabel ( ' r 
(deg/s)'),axis([0 8 -30 30)),grid 

subplot (4,1,4), plot (tld, lr__dr (:, 1) tld, lr_dr (; , 2) , ' - - ' , tld, lr_dr ( :, 3) , ' - ' ) , ylabel ( 'p 

(deg)'),grid 

xlabel('Time (sec)'), ,axis([0 8 -20 20]),pause, 
legend('CMARCClassical','Parameter Est ',0) 
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